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CHAPTER  1 


INTRODUCTION 


1.1  Motivation 

Planetary  gears  are  widely  used  in  the  transmissions  of  helicopters,  automobiles, 
aircraft  engines,  heavy  machinery  and  marine  vehicles.  Figure  1.1  illustrates  a  single- 
stage  planetary  gearset  consisting  of  a  sun  gear,  a  ring  gear,  several  planets,  and  a  carrier. 
Any  of  the  carrier,  ring,  and  sun  can  be  selected  as  the  input  or  output  component,  and  the 
power  is  transmitted  through  multiple  paths  of  the  planet  meshes.  Planetary  gears  have 
substantial  advantages  over  parallel  shaft  drives,  including  compactness,  large  torque-to- 
weight  ratio,  diminished  loads  on  shafts  bearings,  and  reduced  noise  and  vibration  due  to 
the  relatively  smaller  and  stiffer  components. 

Despite  planetary  gears’  advantages,  noise  and  vibration  remain  major  concerns  in 
their  applications.  In  most  helicopters,  planetary  gears  are  used  in  the  last  stage  of  gear 
reduction.  This  planetary  gear  is  mounted  directly  to  the  helicopter  cabin,  so  its  vibration 
is  the  main  source  of  structure-borne  cabin  noise,  which  can  exceed  100  dB  (Krantz, 
1993).  Measurement  of  the  cabin  noise  shows  that  gear  mesh  frequency  and  its  harmonics 
are  the  dominant  acoustic  frequencies  (Figure  1.2).  Extensive  cabin  noise  results  in 
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Figure  1.1  Planetary  gear  diagram  (Lynwander,  1983) 


Frequency,  kH2 

Figure  1.2  Helicopter  cabin  noise  spectra  (Heath  and  Bossier,  1993) 
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operator  fatigue,  communication  difficulty  and  health  risk  from  extended  exposure. 
Mitigating  the  cabin  noise  becomes  a  crucial  requirement  for  effective  operation  of 
military  and  civilian  helicopters.  In  vehicle  automatic  transmissions,  planetary  gear 
vibration  can  be  transferred  through  bearings  and  mounts  to  the  passenger  compartment 
and  results  in  discrete  high  frequency  pitches  in  automotive  interior  noise.  Consumers 
often  perceive  the  gear  noise  as  poor  quality  and  mechanical  problems.  Vibration 
reduction  in  planetary  gears  provides  substantial  benefits:  reduced  noise,  improved 
reliability,  more  efficient  power  transfer,  and  reduced  maintenance  costs.  However,  the 
noise  and  vibration  issues  in  current  planetary  gear  design  rely  mostly  on  empirical 
experience  and  cut-and-try  methods,  rather  than  engineering  understanding.  General 
design  guidelines  are  needed  to  minimize  planetary  gear  vibration. 

Vibration  reduction  of  planetary  gears  requires  thorough  examination  of  the 
structural  dynamics.  A  detailed  review  of  planetary  gear  dynamic  analyses  before  1980s 
is  found  in  August  (1983).  These  analyses  on  planetary  gear  dynamics  include  modeling 
and  free  vibration  investigations  (Cunliffe  et  al.,  1974;  Botman,  1976;  Frater  et  al,  1983; 
Velex  and  Flamand,  1996),  neutralization  of  transmission  error  excitation  (Seager,  1975), 
dynamic  tooth  load  due  to  run-out  errors  (Hidaka,  1979e;  1980),  mesh  stiffness  variation 
(Kasuba  and  August,  1984),  load  sharing  among  planets  (Ma  and  Botman,  1984),  and 
torsional  vibration  and  dynamic  loads  (August  and  Kasuba,  1986).  More  analytical 
studies  were  performed  after  1990s.  Kahraman  derived  a  nonlinear,  time-varying  planar 
model  (1994a)  and  subsequently  extended  it  to  three-dimensions  and  examined  the 
influence  of  planet  phasing  on  dynamic  response  (1994b).  Using  a  three-dimensional 
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model  for  helical  gears,  Kahraman  and  Blankenship  (1994)  investigated  the  load  sharing 
and  mesh  phasing  among  planets.  Kahraman  (1994c)  also  reduced  his  model  to  a  purely 
torsional  one  to  predict  natural  frequencies  and  vibration  modes.  More  recently,  Agashe 
(1998)  and  Parker  et  al.  (2000a)  used  a  finite  element  tool  to  investigate  the  dynamic 
response  and  planet  phasing  issues  in  planetary  gears.  This  special  computational  tool 
naturally  includes  the  time-varying  mesh  stiffness  and  transmission  error  without  ad  hoc 
specification  of  these  factors.  Parker  (2000)  also  rigorously  proved  the  effectiveness  of 
using  planet  phasing  schemes  to  suppress  planetary  gear  vibration. 

The  experimental  studies  on  practical  planetary  gear  vibrations  are  scarce  due  to  the 
difficult  access  to  the  internal  gears.  Chiang  and  Badgley  (1973)  investigated  the  noise 
spectra  generated  from  ring  gear  vibrations  in  the  planetary  reduction  gearbox  of  two 
helicopters  (Boeing- Vertol  CH-47  and  Bell  UH-1D).  Toda  and  Botman  (1979) 
experimentally  showed  that  planetary  gear  vibration  resulting  from  spacing  errors  can  be 
minimized  by  proper  indexing  of  the  planets.  Botman  (1980)  presented  some  typical 
measurement  results  on  the  planetary  gear  of  a  PT6  aircraft  engine.  His  experiments 
showed  some  peculiar  behavior  of  planetary  gear  vibration  regarding  load  sharing, 
response  due  to  gear  errors,  and  dynamic  instability.  Hidaka  and  his  colleagues 
experimentally  studied  the  dynamic  behavior  of  a  Stoeckicht  planetary  gear  and 
published  a  series  of  reports  (Hidaka  et  al,  1976a,b;  1977;  1979a~d).  Their  reports 
studied  some  important  issues  such  as  load  distribution,  effect  of  different  meshing-phase 
among  sun/ring-planet  meshes,  etc.  Velex  et  al.  (1994)  matched  the  natural  frequency 
measurement  of  a  Stoeckicht  epicyclic  train  with  their  finite  element  calculations.  Rakhit 
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(1997)  measured  the  subsynchronous  vibrations  at  the  turbine  bearings  in  a  gas 
turbogenerator  and  proposed  a  design  of  the  epicyclic  gearbox  to  reduce  the  vibration. 
Kahraman  (1999)  developed  a  generalized  model  to  predict  load  sharing  of  planets  under 
quasi-static  conditions  and  validated  the  model  with  experiments. 

According  to  a  comprehensive  literature  research,  less  analytical  investigations 
have  been  done  on  planetary  gear  dynamics  than  those  for  parallel  shaft  gears.  This  is 
largely  due  to  the  modeling  complexity  of  planetary  gears.  Important  complications 
include  multiple  mesh  contacts,  detailed  kinematics,  mesh  stiffness  variation, 
transmission  error  excitation,  contact  loss  nonlinearity,  elastic  ring  gear  vibration  and 
geometric  imperfections.  Most  previous  research  uses  numerical  or  experimental  methods 
to  examine  specific  planetary  gear  systems.  Some  critical  issues  remain  unsolved  and 
require  systematic  analytical  study. 


1.2  Critical  Topics 

The  fundamental  task  of  analytical  planetary  gear  research  is  to  build  a  dynamic 
model.  For  different  analysis  purposes,  there  are  several  modeling  choices  such  as  a 
simple  dynamic  factor  model,  compliance  tooth  model,  torsional  model,  and  geared  rotor 
dynamic  model  (Ozguven  and  Houser,  1988).  According  to  the  source-path-receiver 
relationship  between  the  planetary  gear,  bearing/mounting,  and  the  cabin,  different 
boundaries  can  be  selected  for  building  the  model.  This  study  focuses  on  the 
understanding  of  planetary  gear  dynamic  behavior,  so  a  single  stage  gearset  with  discrete 
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elements  is  the  basic  model  for  investigation.  In  previous  lumped-parameter  models 
(Cunliffe  et  al.,  1974;  Botman  1976;  Kahraman,  1994a,  b,  c),  the  gyroscopic  effects 
caused  by  carrier  rotation  have  not  been  considered.  Because  planetary  gears  have  planets 
mounted  on  the  rotating  carrier,  the  Coriolis  and  centripetal  accelerations  caused  by  the 
carrier  rotation  introduce  gyroscopic  terms  into  the  system  model.  For  high-speed 
applications  such  as  aircraft  engines  (30,000  rpm),  gyroscopic  effects  may  heavily  impact 
the  system  stability  and  behavior.  This  project  will  develop  a  planetary  gear  model 
including  the  gyroscopic  effects,  contact  loss  nonlinearity,  mesh  stiffness  variation,  and 
static  transmission  error  excitation.  Despite  the  use  of  the  term  planetary  gear,  this  model 
is  applicable  for  epicyclic  gears  with  any  configuration  (fixed/floating  sun,  ring,  and 
carrier,  and  non-equally  spaced  planets).  The  model  is  the  fundamental  tool  for  the 
analytical  research. 

The  free  vibration  analysis  calculates  critical  parameters  such  as  natural  frequencies 
and  vibration  modes  that  are  essential  for  almost  all  dynamic  investigations.  Cunliffe  et 
al.  (1974)  numerically  identified  the  planetary  gear  natural  frequencies  and  vibration 
modes  for  a  specific  thirteen  degree-of-freedom  system.  Similar  work  has  been  done  by 
Botman  (1976),  Frater  et  al.  (1983),  and  Kahraman  (1994c)  for  other  example  planetary 
gears,  but  no  systematic  characterization  has  been  obtained.  This  project  reveals  the 
unique  structure  of  natural  frequency  spectra  and  vibration  modes  due  to  the  cyclic 
symmetry  of  planetary  gears  (Lin  and  Parker,  1999a).  All  the  vibration  modes  are 
classified  into  rotational,  translational  and  planet  modes  with  distinctive  properties.  The 
structured  vibration  modes  are  rigorously  characterized  for  general  epicyclic  gears  and 
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validated  by  the  computational  results  from  a  finite  element  model.  These  well-defined 
properties  are  not  valid  when  planets  are  arbitrarily  spaced,  but  still  apply  to  practically 
important  case  of  diametrically  opposed  planets  (Lin  and  Parker,  2000a).  The  free 
vibration  properties  are  very  useful  for  further  analyses  of  planetary  gear  dynamics, 
including  eigensensitivity  to  design  parameters,  natural  frequency  veering,  planet  mesh 
phasing,  and  parametric  instabilities  from  mesh  stiffness  variations. 

Another  key  issue  is  how  design  parameters  affect  the  natural  frequencies  and 
vibration  modes.  During  the  design  process,  model  parameters  are  often  altered  to 
evaluate  alternative  design  choices,  reduce  weight,  and  tune  the  system  frequencies  to 
avoid  resonance.  The  influence  of  design  parameters  on  planetary  gear  natural 
frequencies  was  touched  on  in  a  few  papers,  but  general  conclusions  were  not  presented. 
In  the  plots  of  natural  frequencies  versus  design  parameters,  veering  phenomena  (Leissa, 
1974;  Perkins  and  Mote,  1986)  often  occur  and  obstruct  the  tracing  of  eigenvalue  loci 
under  parameter  changes.  In  the  veering  neighborhood,  where  two  eigenvalue  loci 
approach  each  other  and  then  abruptly  veer  away,  the  veering  vibration  modes  are 
strongly  coupled  and  change  dramatically  (Figure  1.3).  It  is  necessary  to  systematically 
study  natural  frequency  and  vibration  mode  sensitivities  and  their  veering  characters  to 
identify  the  parameters  critical  to  planetary  gear  vibration.  In  addition,  practical  planetary 
gears  may  be  mistuned  by  mesh  stiffness  variation,  manufacturing  imperfections  and 
assembling  errors.  For  some  symmetric  structures,  such  as  turbine  blades,  space 
antennae,  and  multi-span  beams,  small  disorders  may  dramatically  change  the  vibration 
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(b) 

Figure  1.3  (a)  Two  eigenvalue  loci  veer,  (b)  The  associated  vibration  modes  change 
dramatically  through  points  A,  B,  and  C.  The  two  modes  are  strongly  coupled  at  point  B. 
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Figure  1 .4  Stability  chart  of  Mathieu  equation  x  +  (S-2e  cos  2t)x  =  0  (Meirovitch, 

1970).  The  hatched  areas  are  instability  regions. 

modes  and  result  in  mode  localization  (Pierre,  1988;  Cornwell  and  Bendiksen,  1992; 
Happawana  et  al.,  1998).  Vibration  modes  with  dominant  motion  localized  in  one  planet 
lead  to  load  sharing  unbalance,  which  can  severely  undermine  the  power  transfer 
efficiency  and  damage  the  gear  teeth  and  bearings.  This  work  presents  a  thorough 
eigensensitivity  analysis  of  the  natural  frequencies  and  vibration  modes  to  key  model 
parameters  for  both  tuned  (cyclically  symmetric)  and  mistuned  planetary  gears.  The 
parameters  considered  include  mesh/support  stiffnesses,  component  masses,  moments  of 
inertia,  and  operating  speed.  Taking  advantage  of  the  structured  vibration  mode 
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properties,  natural  frequency  sensitivities  are  expressed  in  simple,  closed-form  formulae 
which  relate  the  sensitivities  to  modal  strain  and  kinetic  energy  (Lin  and  Parker,  1999b). 
Well-defined  veering  rules  are  derived  from  these  formulae  and  vibration  mode 
properties  to  predict  veering  and  its  strength.  The  influence  of  design  parameters  is 
examined  through  a  benchmark  example.  The  knowledge  of  natural  frequency  spectra, 
vibration  mode  properties,  eigensensitivity  formulae,  and  special  veering  rules  are 
combined  to  provide  considerable  insight  into  planetary  gear  free  vibration. 

It  is  well  known  that  mesh  stiffness  variation  is  a  major  excitation  source  of  gear 
vibration.  For  spur  gears,  the  time-varying  stiffness  is  caused  by  the  alternating  number 
of  teeth  in  contact.  It  is  a  periodic  function  at  mesh  frequency,  which  is  the  number  of 
tooth  mesh  cycles  per  second.  The  mesh  stiffness  variation  serves  as  a  parametric 
excitation  and  results  in  instability  under  certain  conditions.  Parametric  instabilities  are 
particularly  dangerous  because  they  can  occur  at  excitation  frequencies  well  below  the 
system  natural  frequencies.  The  mesh  stiffness  variation  can  be  further  complicated  by 
the  interaction  of  transmission  error  excitation  (Smith,  1983)  and  contact  loss 
nonlinearity  (Blankenship  and  Kahraman,  1995;  Kahraman  and  Blankenship,  1996, 
1997).  Literature  reviews  of  parametrically  excited  systems  can  be  found  in  the  work  of 
Ibrahim  and  Barr  (1978).  For  a  single  pair  of  gears  excited  by  harmonically  varying 
stiffness,  Bollinger  and  Harker  (1967)  used  the  one  degree-of-freedom  Mathieu  equation 
to  determine  the  instability  regions.  The  instability  conditions  are  often  illustrated  in  the 
plots  of  the  exciting  frequency  versus  the  amplitude  of  varying  stiffness,  as  shown  in 
Figure  1.4.  Benton  and  Seireg  (1978,  1980a)  obtained  the  response  to  mesh  stiffness 
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variation  and  external  excitations  at  integer  multiples  of  the  rotation  speed.  They 
experimentally  verify  the  resonance  region  obtained  from  simulation  and  demonstrated 
the  damage  of  parametric  instability  on  gear  teeth.  Amabili  and  Rivola  (1997)  studied  the 
steady  state  response  and  stability  of  the  single  degree  of  freedom  system  with  time- 
varying  stiffness  and  damping.  Other  researchers  (Benton  and  Seireg,  1981;  Kahraman 
and  Blankenship,  1996;  Nataraj  and  Whitman,  1997;  Nataraj  and  Arakere,  1999)  also 
investigated  gear  parametric  instabilities  using  single  degree-of-freedom  models.  For 
multi-mesh  gear  systems,  it  is  surprising  to  find  little  investigation  on  parametric 
instability  in  the  published  literature.  Although  Benton  and  Seireg  (1980b,  1981)  studied 
a  gear  system  with  two  meshes,  they  uncoupled  the  model  into  two  single  degree-of- 
freedom  equations  with  some  simplifications.  Their  conclusions  on  the  mesh  stiffness 
phasing  effect  contradict  another  investigation  (Tordion  and  Gauvin,  1977)  using  an 
infinite  determinant  analysis  (Bolotin,  1964).  This  conflict  will  be  clarified  using 
perturbation  analysis  (Hsu,  1963,  1965;  Nayfeh  and  Mook,  1979)  and  numerical 
integration  methods.  In  addition,  this  work  extends  parametric  analysis  in  two-stage  gear 
systems  to  planetary  gears  where  parametric  excitations  are  more  complicated  as 
different  contact  ratios  and  phasing  conditions  exist  between  the  sun-planet  and  ring- 
planet  meshes.  August  and  Kasuba  (1986)  and  Velex  and  Flamand  (1996)  numerically 
computed  dynamic  responses  to  mesh  stiffness  variations  for  planetary  gears  with  three 
sequential  phased  planets.  Their  results  showed  the  dramatic  impacts  of  mesh  stiffness 
variation  on  dynamic  response,  tooth  loading,  and  load  sharing  among  planets.  The 
operating  conditions  leading  to  planetary  gear  parametric  instability  have  not  been 
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analytically  investigated.  The  well-defined  vibration  mode  properties  will  be  used  to 
derive  the  operating  conditions  leading  to  planetary  gear  instability.  The  boundaries  are 
expressed  in  simple  forms  and  the  effects  of  contact  ratios  and  mesh  phasing  are 
analytically  investigated.  In  practical  design,  planet  mesh-phasing  schemes  are  often 
applied  to  cancel  or  neutralize  the  excitations  from  transmission  errors  (Seager,  1975; 
Kahraman,  1994a;  Kahraman  and  Blankenship,  1994;  Parker,  2000).  This  study  shows 
that  particular  instabilities  are  eliminated  under  certain  phasing  conditions,  which  can  be 
achieved  by  selection  of  design  parameters  according  to  the  analytical  results.  Dynamic 
response  and  tooth  separations  induced  by  parametric  instability  are  numerically 
examined. 


1.3  Scope  of  Investigation 

The  scope  of  this  project  is  to  advance  the  modeling  and  understanding  of  planetary 
gear  dynamics  and  analytically  examine  certain  critical  factors  affecting  planetary  gear 
noise  and  vibration.  This  research  focuses  on  the  analytical  investigation  of  the  following 
specific  tasks. 

•  Derive  a  lumped-parameter  model  for  spur  planetary  gears,  including  different  planet 
phasing,  gyroscopic  effects,  mesh  stiffness  variation,  and  transmission  error 
excitation.  The  model  is  valid  for  general  epicyclic  gears  with  any  number  of  planets 
and  will  be  the  fundamental  tool  for  further  research. 
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•  Analytically  characterize  the  unique  structure  of  the  natural  frequency  spectra  and 
vibration  modes  of  general  planetary  gears.  The  cases  with  equally  and  arbitrarily 
spaced  planets  (including  diametrically  opposed  planets)  are  considered. 

•  Use  the  vibration  mode  properties  to  obtain  simple,  closed-form  formulae  for  the 
eigensensitivities  to  important  design  parameters.  According  to  these  formulae  and 
characterized  natural  frequency  veering  rules,  the  effects  of  design  parameter  changes 
on  planetary  gear  free  vibration  are  investigated. 

•  Investigate  the  parametric  instabilities  caused  by  multiple  time-varying  mesh 
stiffnesses.  Two-stage  gear  systems  are  examined  first  to  clarify  pervious  conflicts  and 
derive  simple  expressions  of  instability  boundaries.  Then,  the  analytical  method  is 
extended  to  planetary  gear  systems.  The  well-defined  modal  properties  are  used  to 
identify  the  effects  of  contact  ratios  and  mesh  phasing  on  planetary  gear  parametric 
instability. 
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CHAPTER  2 


PLANETARY  GEAR  MODELING 


2.1  Modeling  Considerations 

Lumped-parameter  modeling  is  used  in  this  project  for  dynamic  analysis.  All  gears 
are  considered  as  rigid  bodies  and  component  supports  are  modeled  by  springs.  A  single¬ 
mesh  model  is  shown  in  Figure  2.1  for  a  pair  of  spur  gears.  All  bearings/supports  are 
modeled  as  two  perpendicular  springs.  The  gear  tooth  meshes  are  represented  by  springs 
acting  on  the  line  of  action  with  parallel  viscous  dampers.  The  transmission  error  e(t),  a 
prescribed  displacement  input,  is  included  as  indicated.  The  tooth  separation  nonlinearity 
acts  as  a  step  function  h(t)  where  h=l  when  the  teeth  are  in  contact  and  h-0  when  the 
teeth  lose  contact.  This  model  can  be  extended  to  planetary  gears  with  multiple  meshes. 

In  planetary  gear  modeling,  the  following  factors  must  be  considered. 

1 .  Assumptions.  The  analysis  deals  with  planar  vibration  of  single  stage  planetary  gears. 
Helical  gears  require  three-dimensional  modeling  and  are  not  considered  here. 
Excellent  lubrication  is  assumed  and  tooth  friction  forces  are  neglected.  No  damping 
is  included,  although  a  viscous  damper  could  easily  be  added  in  parallel  with  the 
mesh  and  bearing  springs. 
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Figure  2.1  A  single  mesh  model 

2.  Versatile  configurations.  There  are  many  configurations  for  general  epicyclic  gears. 
By  fixing  one  or  more  of  the  coaxial  components  (carrier,  ring,  and  sun),  various 
configurations  such  as  planetary,  star,  and  differential  gears  can  be  obtained.  In 
typical  designs,  one  of  the  coaxial  members  is  free  to  translate  to  enhance  the  load 
sharing  among  planets.  The  number  of  planets  and  their  positions  also  vary  in 
practical  applications.  The  proposed  model  should  accommodate  these  configurations 
and  be  flexible  for  general  applications. 

3.  System  coordinates.  Several  choices  of  coordinates  may  be  used  in  the  modeling.  The 
frame  can  be  fixed  or  rotating;  the  planet  coordinates  can  be  parallel  to  each  other  or 
use  local  radial  and  tangential  directions.  The  coordinate  selection  does  not  change 
the  physical  properties  of  the  system,  but  may  greatly  affect  the  analysis  difficulty. 
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2.2  System  Equations 

The  planetary  gear  model  used  in  our  analysis  is  shown  in  Figure  2.2.  Each 
component  has  three  degrees  of  freedom:  two  translations  and  one  rotation.  The  model  is 
similar  to  that  used  by  Kahraman  (1994a)  with  two  distinctions:  (1)  the  planet  deflections 
are  described  by  radial  and  tangential  coordinates,  and  (2)  gyroscopic  effects  induced  by 
carrier  rotation  are  modeled.  The  radial  and  tangential  coordinates  more  naturally 
describe  the  vibration  modes.  Gyroscopic  effects  in  high-speed  applications  such  as 
aircraft  engines  may  dramatically  alter  the  dynamic  behavior  from  that  at  lower  speeds. 
The  coordinates  illustrated  in  Figure  2.2  are  used.  The  carrier,  ring  and  sun  translations 
Xh,  yh,  h=c,r,s  and  planet  translations  =  are  measured  with  respect  to  a 

rotating  frame  of  r  reference  i,  j,  k  fixed  to  the  carrier  with  origin  o.  The  Xh,  h=c,r,s  are 
directed  towards  the  equilibrium  position  of  planet  1,  and  £„,r}n  are  the  radial  and 
tangential  deflections  of  the  n-th  planet.  The  basis  i,  j,  k  rotates  with  the  constant  carrier 
angular  speed  Qc-  The  rotational  coordinates  are  uh  =  rhOh,h  =  c,r,s, l,--N ,  where  Oh  is 
the  component  rotation;  is  the  base  circle  radius  for  the  sun,  ring  and  planet,  and  the 
radius  of  the  circle  passing  through  the  planet  centers  for  the  carrier.  Circumferential 
planet  locations  are  specified  by  the  fixed  angles  y/„,  where  %  is  measured  relative  to  the 
rotating  basis  vector  i  so  that  y/i  =  0.  The  details  of  the  model  derivative  are  given  in  Lin 
and  Parker  (1999a).  The  equations  of  motion  are 

Mq  +  QcGq  +  [Kft+Km-^KQ]q  =  T(0  +  F(0  (2.1) 

q  (^c  -  y  c ,  ^  c '  ^  r  *  y  r 7  ^  r‘  7  Tv  ?  ^  s'  ?**’■?  Cn  ?  %  ) 

carrier  ring  sun  planet  1  planet  N 
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where  the  matrix  components  are  given  in  the  Appendix  A.  M  is  the  inertia  matrix  and 
Kb  is  the  bearing  stiffness  matrix.  G  and  Kq  result  from  high-speed  carrier  rotation  and 
have  not  been  included  in  published  models.  To  model  the  time-varying  stiffness 
associated  with  changing  numbers  of  teeth  in  contact  at  each  mesh,  Km  can  be 
decomposed  into  mean  and  time- varying  components.  Tooth  separation  nonlinearity  is 
implicitly  included  in  K m(t).  T (t)  is  the  applied  external  torque  and  F (t)  represents  the 
static  transmission  error  excitation. 

This  comprehensive  model  is  the  fundamental  tool  for  further  research  of  planetary 
gear  dynamics.  It  is  readily  applicable  for  specific  configurations  of  epicyclic  gears. 
Assigning  a  large  value  to  a  transverse  or  torsional  stiffness  restrains  the  translation  or 
rotation  of  a  component;  assigning  a  very  compliant  stiffness  in  a  component  support 
floats  the  component. 


Figure  2.2  Planetary  gear  model 
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CHAPTER  3 


NATURAL  FREQUENCY  AND  VIBRATION  MODE  PROPERTIES 


This  chapter  analytically  investigates  the  natural  frequency  spectra  and  vibration 
modes  of  general  planetary  gears.  The  cyclic  symmetry  of  planetary  gears  leads  to 
highly-structured  free  vibration  characteristics  that  are  identified  herein.  Unique 
properties  of  the  eigensolutions  for  the  linear  time-invariant  case  are  presented  for  an 
example  system.  The  identified  properties  are  then  mathematically  shown  to  be 
characteristics  of  general  planetary  gears. 

3.1  Equally  Spaced  Planets 

The  free  vibration  of  the  linear,  time-invariant  representation  is 

Mq  +  £?tGq  +  (K-fit2KJ<l=®  (3.1) 

where  K=Kb+Km.  The  carrier  speed  is  assumed  to  be  small  and  the  gyroscopic  terms  G 
and  Kn  are  neglected.  The  associated  eigenvalue  problem  of  (3.1)  is 

0)fM</>j  =  K<t>j  (3.2) 

where  C0j  are  natural  frequencies  and  0,  =[pc,p,.,pi,p1,"-piV]r  are  vibration  modes  with 

P/,  =  [xh,yh,uh]T,h  =  c,r,s  for  the  carrier,  ring  and  sun,  and  p„  =  [C„,J1„,un]T  for  the 
planets.  At  this  stage,  the  planets  are  assumed  identical  and  equally  spaced;  all  planet 
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bearing  stiffnesses  are  equal  kp„=kp,  all  sun-planet  mesh  stiffnesses  are  equal  (ksn-ksp); 
and  all  ring-planet  mesh  stiffiiesses  are  equal  kn,=k,T.  With  these  specifications,  planetary 
gears  are  cyclically  symmetric  structures  that  can  be  divided  into  N  identical  sectors.  The 
cyclic  symmetry  of  planetary  gears  leads  to  distinctive  natural  frequency  and  vibration 
mode  properties  that  will  be  demonstrated  analytically. 

We  first  illustrate  the  eigensolution  properties  through  a  numerical  example  with 
the  parameters  in  Appendix  B  case  I.  Typical  vibration  modes  for  N=4  are  shown  in 
Figure  3.1,  where  the  movements  of  the  carrier  and  ring  are  not  shown  in  order  to  clarify 
the  figures.  Some  interesting  conclusions  are  obtained  from  Figure  3.1: 

1.  Six  natural  frequencies  always  have  multiplicity  m=l  for  different  N.  Except  for  the 
zero  natural  frequency,  their  values  increase  as  additional  planets  are  introduced. 
Their  associated  vibration  modes  have  pure  rotation  of  the  carrier,  ring  and  sun 
(Figure  3.1a),  so  these  modes  are  named  rotational  modes. 

2.  Six  natural  frequencies  always  have  multiplicity  m-2  for  various  N.  Some  natural 
frequencies  increase  monotonically  while  others  decrease  monotonically  as  N 
increases.  The  carrier,  ring  and  sun  have  pure  translation  in  the  corresponding 
vibration  modes  (Figure  3.1b),  so  these  modes  are  defined  as  translational  modes. 

3.  Three  natural  frequencies  have  multiplicity  m=N-3  and  exist  only  if  N>3.  Their 
associated  vibration  modes  are  termed  planet  modes  because  the  carrier,  ring  and  sun 
do  not  move;  only  planet  motion  occurs  in  these  modes  (Figure  3.1c,  d,  e).  For  each 
of  these  three  natural  frequencies,  the  corresponding  vibration  modes  span  an  N-3 
dimensional  eigenspace. 
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This  example  shows  there  are  at  most  fifteen  different  natural  frequencies  for  N>3; 
additional  planets  only  change  the  multiplicity  of  the  planet  mode  natural  frequencies. 
When  N<3,  all  natural  frequencies  are  distinct  and  the  vibration  modes  do  not  have 
special  structure  because  of  the  loss  of  cyclic  symmetry.  In  nearly  all  planetary  gear 
designs,  three  or  more  planets  are  used  to  take  advantage  of  the  load  sharing  among 
planets  and  subsequent  discussion  is  restricted  to  this  case.  The  eigensolution  properties 
identified  in  the  example  are  analytically  shown  to  be  true  for  general  planetary  gears 
(Lin  and  Parker,  1999a)  and  summarized  below. 

Planet  Mode.  A  planet  mode  has  the  form 

&  =[0,0,0,  w1p1,.--wjvp1f  (3.3) 

where  w„  are  scalars  (wj=l)  satisfying 

Zw»sin^»=0  Xw„cos^M=0  X>„=°  (3-4) 

Planet  modes  have  the  following  characteristics: 

(i)  The  associated  natural  frequency  has  multiplicity  N-3, 

(ii)  The  translation  and  rotation  of  the  carrier,  ring  and  sun  are  zero. 

(iii)  The  planet  deflections  are  scalar  multiples  of  the  first  (or  any  other  arbitrary) 
planet's  deflection  components. 

Rotational  Mode.  A  rotational  mode  has  the  form 

=[Pc>PrdL>Pl>-,Plf  (3.5) 

Rotational  modes  have  the  following  characteristics: 

(i)  The  associated  natural  frequency  is  distinct, 
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(ii)  The  translation  of  the  carrier,  ring  and  sun  are  zero,  i.e.,  p;,  =  [0,0,M/,]r ,  where  j=c,r,s, 

(iii)  All  planets  have  identical  deflections,  i.e.,  p,  =p2  =•  -^p^ 

Translational  Mode.  A  pair  of  translational  modes  has  the  form 


Pc 

Pc 

P  r 

p, 

P, 

<P,= 

p. 

(cos^pj  Tsin^p,) 

(-  sin  ^p,  +  cos  ^  p, ) 

(cos^p!  4-sin^p,) 

(-sin^p,  +  cos^p,)_ 

Translational  modes  have  the  following  characteristics: 

(i)  The  associated  natural  frequency  has  multiplicity  two. 

(ii)  The  rotation  of  the  carrier,  ring  and  sun  are  zero.  Furthermore,  the  carrier,  ring,  and 
sun  translations  in  the  degenerate  modes  $  and  fa  are  related  by  ph  =  [xh,yh  .0]rand 

P i,  =  [-yh,x,„0]T  ,h=c,r,s, 

(iii)  The  planet  deflections  for  a  pair  of  vibration  modes  are  related  by 


P„ 

cos(//„I  sin  y/lt  I 

Pi 

,P»_ 

-  sin  y/„  I  cos  y/„  I 

.Pi. 

where  p()  and  p„  =  [£„,r}n,un]T  are  the  deflections  of  the  n- th  planet  in  <p,  and  fa .  I  is 
a  3x3  identity  matrix  and  y/n  =  2k  (n - 1 )/ N . 
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(a)  rotational  mode 


(b)  translational  mode 


(e)  rotational  planet  mode 


Figure  3.1  Typical  vibration  modes.  Dashed  lines  are  the  equilibrium  positions 
and  solid  lines  are  the  deflected  positions.  Dots  represent  the  component  centers. 
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3.2  Arbitrarily  Spaced  Planets 

In  general,  much  of  the  above  well-defined  structure  of  the  natural  frequency 
spectra  and  vibration  modes  is  lost  when  the  planets  are  arbitrarily  spaced.  A  notable 
exception  is  the  planet  modes.  Additionally,  for  the  practically  important  case  of 
diametrically  opposed  planets,  the  free  vibration  retains  its  unique  properties. 

Planet  Mode:  As  for  planetary  gears  with  equally  spaced  planets,  systems  with 
arbitrary  planet  spacing  always  have  three  sets  of  planet  modes  of  the  form  (3.3)  with 
multiplicity  N-3.  Only  the  coefficients  w„  obtained  from  (3.4)  are  affected  by  %. 

Rotational  Mode:  In  general,  the  rotational  and  translational  modes  couple  together 
for  arbitrary  planet  spacing  and  no  special  modal  structure  can  be  identified.  For  certain 
planet  spacing,  however,  they  still  have  distinguishing  properties.  A  case  of  particular 
interest  is  that  of  diametrically  opposed  planets,  which  is  common  in  industrial 
applications.  Consider  a  system  with  each  of  N/2  pairs  of  planets  located  along  arbitrarily 
oriented  diameters.  A  pair  of  opposing  planets  have  the  position  relation  %+N/2=W»+n- 111 
this  case, 

sin  VC  =°  Hwn  cos^„  =0  Zu'„  =°  (3-8) 

Accordingly,  systems  satisfying  (3.8)  have  six  rotational  modes  with  property  (3.5).  For 
arbitrarily  distributed  planets  not  satisfying  (3.8),  rotational  modes  do  not  exist. 

Translational  Mode:  While  the  translational  modes  couple  with  the  rotational 
modes  for  truly  arbitrary  planet  spacing,  they  retain  their  structure  for  systems  satisfying 
(3.8).  The  notable  difference  with  equally  spaced  planet  systems  is  that  the  natural 
frequencies  are  no  longer  degenerate  because  the  cyclic  symmetry  is  disturbed.  To  start 
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with,  the  planet  deflection  relations  in  a  translational  mode  are  derived  from  (3.7)  for  any 
three  planets  i,j,  k, 

sin  (Vi  ~  Vj  )P*  +  sin(^.  -  y/k  )p,.  +  sin(^  -  y/t  )py  =  0  (3.9) 

Thus,  the  n-th  planet  deflection  can  be  expressed  as  a  linear  combination  of  pi  and  p2. 
The  component  modal  deflections  for  a  translational  mode  become 

P h=[x„,yh,uh],  h  -  c,r,s,  p„  =[sin(^2-^n)p]  +  sin^„p2]/sin^2  (3.10) 

Therefore,  planetary  gears  with  planet  positions  satisfying  (3.8)  (for  example, 
diametrically  opposed  planets)  have  twelve  distinct  vibration  modes  that  have  the  special 
structure  (3.10)  of  a  translational  mode. 


3.3  Modal  Strain  and  Kinetic  Energy 

The  vibration  modes  can  be  further  characterized  by  modal  strain  and  kinetic 
energy.  The  total  modal  strain  energy  U  and  kinetic  energy  T  are  related  to  the  natural 
frequencies  and  vibration  modes  by 

U  =  =  U,  +  U„  +Ur+Um+U,+  U„  +f1(,U,+Um+UJ  (3.11) 

L  n=\ 

=TI+Tc,+Tr  +  T„+T,+T„l+fJ(T„+T„l)  (3.12) 

Z  n= 1 

where  (i)  £4,  Uhu,  h=  c,r,s  are  the  strain  energies  in  the  translational  and  rotational 
support  springs,  respectively,  of  the  carrier,  ring  and  sun;  (ii)  U„,  Um,  Us„,  n=l,  ...,  N  are 
the  strain  energies  in  the  n-th  planet  bearing,  ring-planet  mesh  and  sun-planet  mesh;  and 
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(iii)  Th,  Thu,  h=c,  r,  s  and  T,„  Tm,  are  the  modal  translational  and  rotational  kinetic 
energies,  respectively,  of  the  carrier,  ring,  sun,  and  planets.  The  detail  of  these  individual 
energies  are 

U  i,  =  (xh  +  yh ),  Ulw  —  —  klmuh ,  h  =  c,r,s 

U„  +(rf.)2],  U.  Um 

Ml  1 

Th=-fco?{xl+yl),  Th=^jCO?ul  h  =  c,r,s  (3.13) 

ni  ] 

T,=-±at(£  +  i£),  Tm=-TjCoful 

Zrp 

where  the  mesh  and  bearing  deflections  8  are  defined  below. 

sun-planet  mesh:  Sm  =yscosy/sll-xssirn//s>,-C„smas-rincosas  +  us+ul,+esll  (3.14) 
ring-planet  mesh:  Sm  =  yr  cos  -xr  sin^ro  ~C„  since,.  +tj„  cos  a,  +ur -u„  +em  (3.15) 
planet  bearing  radial:  8m  =  yc  sin  y/n  +  xc  cosy/,,  -  C,n  (3.16) 

planet  bearing  tangential:  8nt  =  yc  cos y/n  -  xc  sin  \j/n  -  i]n  +  uc  (3.1 7) 

For  the  example  planetary  gear  (Appendix  B,  case  II),  the  modal  strain  and  kinetic  energy 
distributions  are  shown  in  Figure  3.2,  and  the  dominant  motion  and  strain  energy  are 
listed  in  Table  3.1.  According  to  (3.13),  the  dominant  kinetic  energy  occurs  in  the 
component  with  dominant  motion.  Modal  strain  energy  indicates  those  vibration  modes 
most  susceptible  to  parameter  variation  and  identifies  the  most  heavily  loaded  component 
for  response.  In  the  next  chapter,  the  modal  energy  distribution  will  be  used  in  the 
derivation  of  simple  formulae  to  calculate  eigensensitivity  to  key  design  parameters. 
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3.4  Discussions 


Planet  modes  of  multiplicity  N-3  are  remarkably  insensitive  to  planet  location  and 
retain  their  special  properties  for  arbitrary  planet  spacing.  Coupling  between  rotational 
and  translational  modes  occurs  for  arbitrary  planet  spacing,  and  distinct  properties  can 
not  be  identified.  For  systems  satisfying  (3.8),  however,  rotational  and  translational 
modes  have  structured  properties.  This  includes  the  common  case  of  equally  spaced  and 
diametrically  opposed  planet  pairs.  The  foregoing  development  applies  to  general 
epicyclic  gear  configurations.  Configurations  having  fixed  or  floating  components  are 
obtained  by  letting  the  associated  bearing  stiffness  approach  zero  or  infinity  to  obtain  the 
eigensolution  properties.  For  example,  for  the  cases  of  a  fixed  ring,  fixed  sun  or  fixed 
carrier,  a  3(N+2)  degree  of  freedom  system  is  obtained.  The  vibration  modes  in  such 
systems  consist  of  five  rotational,  five  pairs  of  translational  and  three  groups  of  planet 
modes. 

The  free  vibration  properties  have  been  validated  through  a  finite  element 
computation.  A  finite  element  model  (Figure  3.3)  for  this  system  was  built  by  Agashe 
(1998)  and  Parker,  et  al.  (2000)  and  analyzed  with  CALYX  (Vijayakar,  1991). 
Computational  modal  analyses  were  performed  by  applying  impulse  inputs  and  obtaining 
the  frequency  response  functions.  The  identified  natural  frequencies  match  the  analytical 
results  (Table  3.2)  within  3.3  percent  difference. 
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(a)  Normalized  modal  strain  energy  at  the  sun/planet  meshes  (1-4),  ring/planet  meshes 
(5-8),  planet  bearings  (9-12),  carrier  and  sun  bearings  (13,  14). 
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(b)  Normalized  modal  kinetic  energy  at  the  carrier,  sun,  planet  translations  (1,2, 3-6),  and 
carrier,  sun,  planet  rotations  (7,8,9-12). 


Figure  3.2  Modal  energy  distributions  of  the  system  in  Appendix  B,  case  II. 
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Type 

Dominant  motion 

Dominant  strain  energy 

1 

0 

R 

Rigid  body  mode 

2,3 

825 

T 

All  components 

Carrier,  sun  bearings 

4 

1661 

R 

All  components 

Planet  bearings 

5 

1808 

P 

Planet  radial  translation 

Planet  bearings 

| 

1834 

T 

Sun,  planet  translation 

Sun,  planet  bearings 

8 

1985 

R 

Sun  rotation,  planet  translation 

Planet  bearings 

9,10 

2326 

T 

Sun,  planet  translation 

Sun,  planet  bearings 

11 

5964 

P 

Planet  tangential  translation 

All  meshes 

12,13 

6429 

T 

Planets 

Ring-planet  meshes 

14 

6451 

R 

Planets 

Ring-planet  meshes 

15 

6982 

P 

Planet  rotation 

All  meshes 

16,17 

10430 

T 

Sun  translation 

Sun-planet  meshes 

18 

13068 

R 

Sun  rotation 

Sun-planet  meshes 

Table  3.1  Dominant  motion  and  modal  strain  energy  in  the  vibration  modes  of  the 
system  in  Appendix  B.  R:  rotational  mode,  T:  translational  mode,  P:  planet  mode. 


Mode  Type 

Trans 

Rot 

Plan 

Trans 

Rot 

Trans 

FEM  (Hz) 

778 

1144 

1729 

1676 

Analytical  (Hz) 

769 

1156 

1609 

1710 

1781 

2175 

Difference  (%) 

1.2 

-1.0 

2.3 

-2.0 

-3.3 

-3.0 

Table  3.2  Comparison  of  analytical  and  FEM  natural  frequency  analyses. 
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Figure  3.3  Finite  element  model  of  the  planetary  gear 


The  well-defined  eigensolution  properties  are  useful  for  subsequent  research 
because  almost  all  analytical  investigation  of  planetary  gear  vibration  phenomena 
ultimately  require  the  natural  frequencies  and  vibration  modes.  This  includes,  for 
example,  critical  dynamic  behaviors  such  as  forced  response  to  static  transmission  error, 
use  of  planet  phasing  to  eliminate  excitation  of  particular  modes,  parametric  instability 
from  time-varying  mesh  stiffness,  contact  loss  nonlinearities,  natural  frequency  and 
vibration  mode  sensitivity  to  key  design  parameters  and  gyroscopic  effects.  The  specific 
properties  characterized  in  this  paper  theoretically  explain  the  selective  participation  of 
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the  different  classes  of  vibration  modes  in  planetary  gear  dynamic  response  (Parker  et  al., 
2000a;  Parker,  2000).  Additionally,  the  structured  vibration  mode  properties  lead  to 
simple,  exact  formulae  to  calculate  natural  frequency  and  vibration  mode  sensitivity  to 
parameter  changes  (Lin  and  Parker,  1999b).  Using  these  properties,  well-defined  veering 
rules  of  natural  frequency  loci  are  analytically  derived  (Lin  and  Parker,  2000b).  These 
special  veering  rules  help  to  trace  the  evolution  of  the  loci  for  changes  in  the  design 
parameters  and  identify  the  veering  zones  where  vibration  modes  undergo  dramatic 
changes.  Apart  from  analytical  applications,  the  identified  natural  frequency  structure 
provides  important  information  for  tuning  the  system  frequencies  to  avoid  resonance. 
There  are  at  most  15  different  natural  frequencies  in  general  planetary  gears  with  N 
planets,  so  there  are  only  15  potential  resonant  frequencies. 
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CHAPTER  4 


EIGENSENSITIVITY  TO  DESIGN  PARAMETERS 


During  the  design  process,  system  parameters  are  varied  to  evaluate  alternative 
design  choices,  avoid  resonances,  optimize  load  distribution,  and  reduce  weight.  It  is 
important  to  characterize  the  effects  of  parameter  variations  on  the  natural  frequencies 
and  vibration  modes  for  effective  vibration  tuning.  In  planetary  gear  dynamic  models 
(Figure  2.2),  the  key  design  parameters  include  the  mesh  stiffnesses,  support/bearing 
stiffnesses,  component  masses,  and  moments  of  inertia.  The  influence  of  some  design 
parameters  on  planetary  gear  natural  frequencies  was  touched  on  in  a  few  papers.  Botman 
(1976)  and  Cunliffe  et  al.  (1974)  both  presented  plots  of  natural  frequencies  versus  planet 
bearing  stiffness.  Kahraman  (1994c)  showed  the  effects  of  mesh/bearing  stiffnesses  on 
the  natural  frequencies  in  his  torsional  model  of  planetary  gears.  Saada  and  Velex  (1995) 
discussed  the  influence  of  ring  support  stiffness  on  free  vibration.  These  analyses  were 
based  on  parametric  studies  of  example  planetary  gears  and  assume  the  planetary  gears  to 
be  cyclically  symmetric  (tuned)  systems.  Eigensensitivity  analysis  for  mistimed  systems 
is  necessary  to  identify  the  critical  modes  that  are  susceptible  to  irregularity.  Frater  et  al. 
(1983)  studied  the  vibration  modes  with  one  imbalanced  mesh  stiffness,  but  general 
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conclusions  were  not  obtained.  This  project  analytically  investigates  the  natural 
frequency  and  vibration  mode  sensitivity  to  most  system  parameters  such  as 
mesh/bearing  stiffnesses,  component  masses  and  moments  of  inertia.  Simple,  closed-form 
expressions  are  obtained  to  calculate  eigensensitivities  for  both  tuned  and  mistuned 
system.  In  addition,  eigenvalue  veering  phenomena  are  investigated  to  identify  dramatic 
changes  of  natural  frequencies  and  strong  coupling  of  vibration  modes.  Design  guidelines 
are  summarized  from  eigensensitivity  and  veering  analyses  to  predict  influences  of 
system  parameters  on  planetary  gear  free  vibration. 


4.1  Calculation  of  Eigensensitivity 

The  eigensensitivity  analysis  calculates  natural  frequency  and  vibration  mode 
derivatives  with  respect  to  stiffnesses,  masses,  moments  of  inertia  and  the  carrier  rotation 
speed  Qc-  Eigensensitivity  to  stiffness  and  inertia  design  parameters  are  examined  in  the 
absence  of  gyroscopic  effects  Qc.  Gyroscopic  effects  are  important  in  high-speed 
applications  such  as  aircraft  engines,  and  eigensensitivity  with  respect  to  Qc  is  studied 
separately  in  section  4.4.  The  eigenvalue  problem  for  the  study  is  in  form  (3.2),  i.e., 

(K-A,MM=0  (4.1) 

where  Af=C0j2.  The  eigensensitivity  for  problems  in  the  form  (4.1)  has  been  thoroughly 
investigated  (Courant  and  Hilbert,  1953;  Adelman  and  Haftka,  1986;  Friswell,  1996)  and 
the  necessary  results  are  introduced  below.  The  unique  modal  properties  of  planetary 
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gears  are  then  invoked  to  reduce  these  general  results  to  simple  expressions  specific  to 
planetary  gears. 

Let  (  )'  and  (  )"  denote  the  first  and  second  derivatives  with  respect  to  a  model 
parameter  (i.e.  mesh/bearing  stiffness,  component  mass  or  moment  of  inertia).  For 
simplicity,  the  eigenvalue  derivatives  X{  and  X"  are  calculated;  the  relations  (dj '=Xj'/(2cQj) 
and  (Oi  "=(2AiAj  '- A*  z)/(4 o^)  yield  the  natural  frequency  sensitivities.  For  a  distinct 
eigenvalue,  the  eigensensitivities  are  (Fox  and  Kapoor,  1968;  Rogers,  1970) 


L  AT  , 


1  _  vL  d>J (K  '  -  A,.M  ')<bi 

</>'=  +  Z  — - — — 

2.  k-i  A-  Ai, 


(4.2) 

(4.3) 


2 <j)f  (K  '  -  A,M  ')</><  +  <t>  f  (K  "  -  A,M  "  -  2  A/M  ')<pt 


(4.4) 


For  the  case  of  degenerate  eigenvalues,  consider  a  system  having  a  group  of 
eigenvalues  A]=...=Xm  with  multiplicity  m.  The  first-order  eigenvalue  derivatives  A, 'are 
the  eigenvalues  of 

Da,  =  /1/a.,  D  =  rT(K'-^M')r  (4.5) 


where  T=[% is  an  arbitrary  set  of  independent  eigenvectors  associated  with  this 
degenerate  eigenvalue  and  is  normalized  such  that  T  MT  =  Imxm .  For  the  case  when  all 
A, 'obtained  from  (4.5)  are  distinct,  the  ai  are  uniquely  obtained  with  the  normalization 
aiTaj=l.  This  procedure  determines  the  set  of  independent  eigenvectors  $=Tai  that  admit 
continuous  change  of  the  eigenvectors  as  the  degenerate  eigenvalues  split  into  distinct 
ones  when  a  parameter  is  varied. 

The  eigenvector  derivatives  for  distinct  X{  are  expressed  as  (Friswell,  1996) 
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#  =  v,-+4>c„ 


(4.6) 


where 


v, .  £  * <pt 

k=m+ 1  \  ~~  A 

(4.7) 

2<p]  (K'- A,M')v,  +(/)]  (K'-A,M'-2A'M')0, 

a;- a; 

j*i 

(4.8) 

c  =  --<DrM'0 
"  2 

(4.9) 

The  second  derivatives  are 

A"=2tf(K'~  A,M>,.  +  ^(K"-XIM"-2X'M%,  (4.10) 

For  the  case  when  all  A, 'obtained  from  (4.5)  are  degenerate,  the  a,  are  not  unique 
and  hence  i=l,...,  m  are  arbitrary  in  the  eigenspace.  The  eigenvector  derivatives  can 
not  be  determined  when  these  degenerate  modes  do  not  separate.  However,  X, "  can  be 
obtained  from  the  eigenvalues  of  (Friswell,  1 996) 

E  =  2<Dr(K'-A;M/)V  +  Or(K'-A,M'-2T;M,)^>  (4.11) 

where  V=[vi,...,vm]  is  determined  by  (4.7).  X, "are  not  affected  by  the  selection  of  O. 

The  foregoing  development  is  used  subsequently  to  derive  general,  closed-form 
eigensensitivity  relations  for  Xj (f)'  and  A, "  for  planetary  gears.  These  expressions  yield 
eigensolution  approximations  according  to 

d?i  d(b 

A  =  A  +  ^]— |/>  =  p0(p- p0),  <f>  =  <[>+'Yj~^r\p  =  p0(P~  Po)  (4.12) 

p  dp  p  dp 

where  p  represents  any  system  parameter  with  nominal  value  p§  and  multiple  parameter 
perturbations  are  permitted.  Eigensolutions  X,  <p  are  for  a  nominal  set  of  model 
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parameters  referred  to  as  the  unperturbed  system,  and  the  derivatives  are  evaluated  for 


this  unperturbed  system.  Eigensolutions  A,  (p  are  for  the  perturbed  system  with  varied 
parameters.  The  unperturbed  system  is  taken  to  be  tuned  in  this  study.  Note  that  this  does 
not  meaningfully  restrict  the  results  because  parameter  variations  leading  to  both  tuned 
and  mistuned  perturbed  systems  are  examined. 


4.2  Eigensensitivity  to  Mesh  and  Support  Stiffnesses 

The  stiffnesses  under  consideration  (Figure  2.2)  include  mesh  stiffnesses  km,  ks„, 
transverse  support  stiffnesses  kc,  kr,  ks,  k„,  and  rotational  support  stiffnesses  kcu,  kru,  ksu. 
The  natural  frequency  sensitivity  to  a  certain  stiffness  is  found  to  be  uniquely  associated 
with  the  modal  strain  energy  occurring  in  that  spring.  To  demonstrate  the  procedure,  let 
the  sun-planet  mesh  stiffness  ks„  be  the  varied  parameter. 


4.2.1  Tuned  System 

Considered  the  case  where  all  sun-planet  mesh  stiffnesses  ks„=ksp  are  altered 
equally  so  the  perturbed  system  remains  tuned.  For  rotational  modes,  the 
eigensensitivities  are  obtained  from  (4.2)-(4.4) 


d\ 

& sp 


&sp 


N 


=n 


8‘SL 


L-j  i 

*= 1  «= 1  A 

kid 


'fa 


(4.13) 

(4.14) 
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(4.15) 


The  rotational  mode  property  (3.5)  dictates  that  all  sun-planet  mesh  deformations  are 
equal,  i.e.,  S'm  =  ,  so  (4.13)  becomes 


d\_ 

dksp 


(4.16) 


In  (4.14),  fa '  is  expressed  as  a  modal  expansion  of  eigenvectors,  and  the  contribution  of 
each  eigenvector  is  readily  obtained  from  the  coefficients  of  fa.  When  two  eigenvalues  A,- 
and  A*  are  nearly  equal,  the  influence  of  fa  on  fa' is  dominant  because  the  denominator  in 
its  coefficient  is  small.  In  such  cases,  the  second  derivative  fa" is  also  large,  and  the 
natural  frequency  changes  rapidly  with  ksp. 

The  translational  mode  eigenvalues  Au  do  not  separate  because  the  perturbed 

system  remains  tuned.  Thus,  the  matrix  D  in  (4.5)  has  degenerate  eigenvalues  A/(  A?( 
Accordingly,  the  unperturbed  eigenvectors  fa  2  can  not  be  uniquely  determined  from  the 
procedure  associated  with  (4.5),  and  fa>\  are  an  arbitrary  pair  of  translational  modes 

of  the  unperturbed  system.  From  (4.5),  X\,  X2  are  the  eigenvalues  of 

s'j;n  1 

D  =  0rK'0  =  Y  "  ,  “  “  (4.17) 

ZmmJ  cl  c2  /  C2  \2  '  ^ 

"=l|AAn  \Psn)  J 
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Use  of  the  translational  mode  property  (3.6)  yields  =T , 

Tl.SlSj  =0.  Thus,  the  eigenvalues  of  D  (i.e.,  Aj '  A21)  are  degenerate  and  have  the 


form  (4.13)  for  i=l,2.  From  (4.1 1)  and  the  translational  mode  properties,  Aj  ”,  A2”  are 


d2\ 


=yy 


2(%A,)2 


La  L-j  2—7  5 

aKsp  k-3  n=l  \ 


Uj  =  1,2 


(4.18) 


Planet  modes  are  also  degenerate  and  the  procedure  is  similar  to  that  for 
translational  modes.  For  planet  modes  ...,(fh\,  the  elements  of  the  matrix  D  are 

Dtj  -  f°r  Applying  planet  mode  property  (3.3)  to  calculate 


results  in  £);/=...=  Dmm  and  Dp=0,  irf.  It  follows  that  all  V  of  a  group  of  planet  modes 
are  equal  and  can  also  be  expressed  as  (4.13)  for  ij=l,...,m.  In  the  same  way,  all  At", 
ij=l,  ...,m  are  equal  and  of  the  form  (4.18). 

Equation  (4.13),  which  is  valid  for  all  three  types  of  vibration  modes,  can  be  related 
to  the  modal  strain  energy  U  in  <)>,•.  With  the  definition  of  strain  energy  Us„  in  (3.13), 
(4.13)  becomes 


d\ 

dksP 


N 

K'sp  H=1 


sn  5 


dC0t  1  A 

OKsp  Wi^sp  n= 1 


(4.19) 


Equation  (4.19)  allows  one  to  obtain  the  natural  frequency  sensitivity  to  sun-planet  mesh 
stiffness  by  inspection  of  the  modal  strain  energy  distribution. 

As  an  example,  consider  a  planetary  gear  used  in  the  transmission  of  a  U.S.  Army 
OH-58  helicopter.  The  nominal  model  parameters  are  listed  in  Appendix  B,  case  II.  The 
natural  frequencies  from  (4.1)  are  shown  in  Figure  4.1a  for  a  range  ksp.  The  strain 
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energies  of  each  spring  are  calculated  according  to  (3.11)  and  their  distribution  in  mode 
16  (a  translational  mode)  is  shown  in  Figure  4.1b,c  for  two  cases:  ksp=70  Nl\x.m  and 
ksp=500  N!\x.m.  The  associated  vibration  mode  (pi  6  is  also  shown  for  these  two  cases. 
Little  strain  energy  is  stored  in  the  sun-planet  meshes  Us„  for  case  I,  while  substantial 
strain  energy  results  in  case  EL  Consequently,  (Ou,  is  more  sensitive  to  kSJ)  in  case  II  than 
in  case  I.  This  conclusion  is  consistent  with  the  larger  slope  of  the  (Oi6  locus  for  case  El  in 
Figure  4.1a.  In  fact,  natural  frequency  sensitivities  to  all  stiffnesses  can  be  obtained 
quantitatively  directly  from  the  strain  energy  distribution  using  (4.19)  and  analogous 
relations  (4.20)  and  (4.21)  below. 


2  2  2 

^rx"+y"-Ku"- 


3-  =  ui=7~U> 


h  =  c,r,s 


OKh  n=l  Kp  »= 1 


^j(^7i)  i  n 

^  rp  n= 1  *rp  «=1 


(4.20) 

(4.21) 


These  relations  apply  for  all  three  types  of  vibration  modes.  Expressions  for  (p, ' and  Xj" 
for  all  of  the  stiffiiess  parameters  are  collected  in  Lin  and  Parker  (1999b).  Recalling  the 
special  properties  of  vibration  modes,  (4. 1 9)-(4.2 1)  imply  that 

1 .  Rotational  modes  are  independent  of  the  transverse  support  stiffnesses  of  the  carrier, 
ring,  and  sun  because  these  components  have  no  deformation  of,  and  hence  no  modal 
strain  energy  in,  their  transverse  support  springs.  Thus,  use  of  a  "floating"  sun,  ring, 
or  carrier  i.e.,  low  stiffness  support)  has  no  impact  on  rotational  modes. 

2.  Translational  modes  are  similarly  independent  of  the  rotational  support  stiffnesses  of 
the  carrier,  ring,  and  sun. 

3.  Planet  modes  are  insensitive  to  all  carrier,  ring,  and  sun  support  stiffnesses. 
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Usn  Um  Un  Uc  Us  Usn  Um  Un  Uc  Us 


Figure  4.1  (a)  natural  frequency  versus  the  sun-planet  mesh  stiffness  ksp.  (b),  (c)  Mode  16 
strain  energy  distribution  in  case  I  and  case  n.  All  U  are  defined  in  (3. 1 1). 
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4.2.2  Mistuned  System 


In  practical  planetary  gears,  mistuning  may  be  caused  by  differing  mesh  stiffnesses 
between  planets  due  to  differing  numbers  of  teeth  in  contact,  manufacturing  variations, 
and  assembly  errors.  To  study  the  effects  of  mistuning  on  eigensolutions,  we  examine  the 
sensitivity  to  parameter  variations  that  differ  between  the  planets.  Consider  an  example 
with  only  the  first  sun-planet  mesh  stiffness  ksi  varying  from  the  nominal  (unperturbed) 
value  ksp.  The  derivatives  of  the  mass  and  stiffness  matrices  with  respect  to  ksi  are 

m'=m"=k"=o,  k'  =  ac  /  dkA . 

The  eigensensitivities  of  the  rotational  modes  are  obtained  from  (4.2)-(4.4) 

(4.22) 

(4.23) 


3-=(#*=±„ 
<*,.  (,l)  ■' 

dty,  y1  4i4'i  . 

dksl~hA-\*k 

k*i 

)2  1  L 


=  2 


-oM.)2 


(4.24) 


dk;,  ft?; 

k*i 

Equation  (4.22)  relates  X,'  to  the  modal  strain  energy  in  the  first  sun-planet  mesh. 
Equations  (4.22)-(4.24)  are  similar  to  (4.13)-(4.15)  without  the  summation  over  n 
because  the  varying  parameter  is  located  only  at  the  first  sun-planet  mesh. 

For  translational  modes,  the  eigensensitivities  are  (Lin  and  Parker,  1999b) 


^1  dX^  .  „2\2 
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The  behavior  of  (pi  and  (jh  is  shown  in  Figure  4.2  for  the  example  of  Appendix  B,  case  II. 
A  pair  of  translational  mode  natural  frequencies  separate  as  a  disorder  £=ksi/ksp-l  is 
introduced.  The  modal  strain  energy  distributions  in  the  four  sun-planet  meshes  are 
shown  for  £  =  0,  -0.1.  (pi  is  sensitive  to  ksi  because  of  the  high  strain  energy  in  the  first 
sun-planet  mesh,  fa  has  no  strain  energy  in  the  first  sun-planet  mesh  and  is  independent 
of  ksj.  The  linear  (  ^  +  £kspAJ  )  and  quadratic  ( A,  +  £kxpX^  +  £2kspA^' )  approximations  of 

the  loci  are  shown  in  Figure  4.2  and  agree  well  with  the  exact  loci.  These  two  loci 
intersect  exactly  at  e  -  0  when  there  is  only  one  disorder  in  the  perturbed  system.  If  one 
more  disorder  £2~kS2/ksp-l =  0.1  is  added  at  the  second  mesh  (Figure  4.3),  the  two  loci 
suddenly  change  direction  and  veer  away.  For  an  initially  tuned  (cyclically  symmetric) 
system,  two  independent  varying  parameters  (e.g.,  ksi  and  ksj)  are  necessary  to  break  the 
symmetry  of  both  (pi  and  (p2  and  cause  frequency  loci  veering  (Happawana  et  al.,  1998). 
The  regular  perturbation  does  not  give  a  good  approximation  in  the  veering  zone  where 
two  loci  are  close  to  each  other.  Singular  or  improved  perturbation  methods  are  needed 
(Pierre  and  Murthy,  1992;  Wu,  1993;  Chen  et  al.,  1995;  Lin  and  Lim,  1997).  If  the 
number  of  planets  N=  4  or  5,  the  planet  modes  have  multiplicity  m-1,2  and  their 
eigensensitivities  are  obtained  from  (4.22)-(4.24)  or  (4.25)-(4.27).  When  N>5, 
eigensolutions  of  matrix  D  in  (4.5)  are  difficult  to  achieve  in  closed-form,  but  can  be 
obtained  numerically. 
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Figure  4.2  Influence  of  the  disorder  e  on  the  natural  frequencies.  Linear  (...)  and 
quadratic  ( — )  approximations  agree  well  with  the  exact  loci  ( — ■). 


xIO4 


Figure  4.3  Influence  of  two  disorders  on  the  natural  frequencies.  Another  disorder 
£2=ks2/ksp- 1=0.1  at  the  second  mesh  is  added  to  the  system  shown  in  Figure  4.2. 
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4.3  Eigensensitivity  to  Gear  Mass  and  Inertia 


The  parameters  of  interest  consist  of  masses  (mc,  mr,  ms,  mp)  and  moments  of  inertia 
(Ic,Ir,  Is,  Ip)  for  the  carrier,  ring,  sun,  and  planets.  When  the  perturbed  system  remains 
tuned,  the  eigenvalue  derivatives  for  the  three  types  of  modes  are 
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(4.29) 


where  Th,  Thu,  h=c,  r,  s  and  T„,  T„u  are  the  modal  kinetic  energies  defined  in  (3.12). 
Expressions  for  </>'  and  for  all  of  the  mass  and  inertia  parameters  are  collected  in  Lin 
and  Parker  (1999b).  Figure  4.4a  shows  an  example  plot  of  the  natural  frequencies  versus 
the  sun  moment  of  inertia  Is.  Most  natural  frequencies  are  insensitive  to  changes  in  Is.  The 
kinetic  energy  distribution  and  vibration  modes  of  mode  18  (a  rotational  mode)  are 
shown  in  Figure  4.4b, c  for  cases  I  and  II.  The  sun  has  more  rotational  kinetic  energy  Tsu 
in  case  I  than  in  case  II,  so  (0]s  locus  has  larger  slope  in  case  I.  Equations  (4.28)  and 
(4.29)  allow  quantitative  calculation  of  natural  frequency  sensitivity  to  all  masses  and 
moments  of  inertia  directly  from  the  modal  kinetic  energy  distributions.  Considering  the 
properties  of  vibration  modes,  some  conclusions  are  immediate  from  (4.28)  and  (4.29): 

1.  Rotational  modes  are  independent  of  the  masses  of  the  carrier,  ring,  and  sun  because 
such  modes  have  no  translations  of  these  components. 

2.  Translational  modes  are  similarly  independent  of  the  moments  of  inertia  of  the 
carrier,  ring,  and  sun. 

3.  Planet  modes  are  independent  of  both  masses  and  inertias  of  the  carrier,  ring,  and  sun. 
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(b)  case  I  (c)  case  II 
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Figure  4.4  (a)  Natural  frequency  versus  the  sun  moment  of  inertia  Is.  (b),  (c)  Mode  18 
kinetic  energy  distribution  in  case  I  and  II.  The  T are  defined  in  (3.12). 
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4.4  Eigensensitivity  to  Operating  Speed 


In  high-speed  applications  ( e.g .  aircraft  engines),  gyroscopic  effects  may 
significantly  alter  the  system  stability  and  dynamic  behavior.  Eigenvalue  derivatives 
evaluated  for  Q 'c=0  are  calculated  to  assess  the  influence  of  operating  speed  on  the 
natural  frequency  spectrum.  The  gyroscopic  eigenvalue  problem  of  (3.1)  is  obtained  from 
the  separable  solution  q  =  §leia),t 

\-co] M  +  jmfifi  +  (K  - Q2cKnM  =  0  .  (4.30) 

For  practical  operating  speeds  i2c= 0  (i.e.,  subcritical),  the  eigenvalues  remain  purely 
imaginary.  Suppose  a  zero  speed  natural  frequency  0),  has  multiplicity  m  and  the 
arbitrarily  chosen  independent  eigenvectors  are  P with  normalization 
rrMr  =  lmxm .  While  eigenvectors  for  12^0  are  complex,  the  y  are  real.  Differentiation  of 
(4.30)  with  respect  to  Qc  and  evaluation  at  Oc=0  yield 

(K  -  o)f  M)<t>;  =  (2coi(o^A-j(oiG)Tai  =  f  (4.31) 

where  =  Ta, .  Applying  solvability  and  normalization  conditions  results  in  an  mxm 
Hermitian  eigenvalue  problem 

Da,.  =«>,.,  D  =  yrrGr  /  2  (4.32) 

The  natural  frequency  sensitivities  6)'  are  obtained  from  the  eigenvalues  of  (4.32)  for  the 
three  classes  of  vibration  modes.  Rotational  mode  natural  frequencies  are  distinct  and 
(4.32)  becomes  a  scalar  equation.  Hence  eo'=  0  because  yjGyt  =  0  for  real  y*  and  skew- 
symmetric  G.  For  translational  modes  yi  and  72,  D  and  its  eigenvalues  are 
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Figure  4.5  Natural  frequency  versus  the  carrier  rotation  speed 
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(4.33) 


For  a  group  of  planet  modes  yj,...,ym,  the  properties  (3.3)  guarantee  DtJ  =  yj  Gy j  =  0  for 

irf\  Du—0,  i=l,...,  m  because  of  the  skew-symmetry  of  G.  Thus,  D=0  and  all  planet 
mode  natural  frequency  sensitivities  vanish,  i.e.,  co'=  0. 

Equation  (4.33)  can  be  used  to  approximate  the  frequency  loci  5)l  =0),  +  QcO)' . 
The  result  co'=  0  for  rotational  and  planet  modes  at  Qc  =  0  indicates  the  natural 
frequencies  of  these  modes  are  scarcely  affected  by  operating  speed.  Figure  4.5  shows  the 
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first  10  frequency  loci  versus  Qc  for  the  gear  system  in  Appendix  B  case  II.  The 
rotational  mode  (0)4,  cog )  and  planet  mode  (0J5)  loci  are  nearly  flat  lines  and  not  sensitive 
to  operating  speed.  Translational  mode  frequencies  (obj,  Ob, 7,  Ob, 10)  split  as  Qc  is 
increased  from  zero.  In  this  example,  Gb  and  0)3  at  Qc=600  rad/s  deviate  about  10  percent 
from  the  zero  speed  value.  Typical  helicopter  carrier  speeds  are  less  than  100  rad/s.  For 
applications  with  high  speed  (e.g.  turbofan  and  turboprop  engine  systems),  heavy 
component  masses,  and  compliant  stiffnesses,  the  gyroscopic  effects  can  be  more 
significant.  If  a  natural  frequency  locus  has  large  slope  and  decreases  to  zero  in  the  range 
of  operating  speed,  the  stability  and  system  behavior  are  dramatically  impacted. 


4.5  Natural  Frequency  Veering 

In  the  plots  of  natural  frequencies  versus  design  parameters,  eigenvalue  veering 
occurs,  where  two  eigenvalue  loci  approach  each  other  as  a  parameter  is  varied  but  then 
abruptly  veer  away  like  two  similar  charges  repelling  (point  B  in  Figure  4.6a).  The 
phenomenon  has  been  studied  extensively  (Leissa,  1974;  Perkins  and  Mote,  1986;  Pierre, 
1988;  Chen  and  Ginsberg,  1992).  The  vibration  modes  of  the  veering  eigenvalues  are 
strongly  coupled  and  undergo  dramatic  changes  in  the  veering  neighborhood.  In  the  case 
of  especially  sharp  veering,  it  is  sometimes  difficult  to  distinguish  between  intersection 
and  veering  just  by  observing  eigenvalue  plots.  When  multiple  curves  veer  or  intersect 
close  together  (Figure  4.6a),  strong  modal  coupling,  and  the  associated  operating 
condition  response  changes  that  occur,  are  not  identifiable  from  frequency  loci  plots.  The 
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objective  of  this  work  is  to  analytically  characterize  the  rules  of  eigenvalue  veering  in 
planetary  gear  vibration.  Simple  rules  emerge  to  predict  if  two  eigenvalues  veer  or  cross. 
The  veering  sharpness  is  also  calculated. 


4.5.1  Veering/Crossing  Criterion 

A  method  for  detecting  eigenvalue  veering/crossing  in  general  dynamic  systems  is 
developed  by  Perkins  and  Mote  (1986).  When  two  eigenvalue  loci  veer  away,  their  loci 
curvatures  indicate  the  abruptness  of  curve  direction  changes.  Perkins  and  Mote 
estimated  the  loci  curvature  in  the  veering  neighborhood  using  coupling  factor 


=  2[^(KHMUf 

Zr  a -  a. 


__  2[0/(K'-TrM/)<U2 
Zs  A,  -  A. 


(4.34) 


The  coupling  factors  %s  approximate  the  local  curvatures  and  the  coupling  strength. 
Figure  4.6b  shows  X=Xu  =  ~Xis  versus  the  varying  parameter  ksp  for  the  veering  loci  6)u 
and  (Oi8  in  Figure  4.6a.  Notice  the  sharply  changing  vibration  modes  indicated  in  Figure 
4.6b.  The  two  veering  loci  exchange  mode  shapes  from  point  A  to  C,  even  though  the  loci 
do  not  intersect.  The  modes  are  strongly  coupled  at  B  and  do  not  look  like  either  of  the 
veering  modes  just  outside  the  veering  zone.  When  a  parameter  is  adjusted  in  the  veering 
zone,  the  drastic  changes  in  the  vibration  modes  can  greatly  impact  the  operating 
condition  dynamic  response,  tooth  loads,  load  sharing,  and  bearing  forces  and  possibly 
lead  to  mode  localization.  The  degree  to  which  individual  modes  are  excited  by  dynamic 
mesh  forces  (i.e.,  the  modal  forces)  also  changes  dramatically  as  veering  alters  the 
modes.  If  the  coupling  factors  are  all  zero,  X  and  As  loci  cross;  otherwise,  veering  occurs. 
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Figure  4.6  (a)  Natural  frequencies  versus  the  sun-planet  mesh  stiffness.  Natural 
frequencies  are  numbered  under  the  nominal  conditions  (dashed  line)  in  Table  C.l.  Three 
pairs  of  veering  are  loci  14  and  18,  (12,13)  and  (16,17),  and  11  and  15.  (b)  Coupling 
factor  x  of  (Ou  and  (Ois-  These  two  loci  exchange  mode  shapes  from  points  to  C. 
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4.5.2  Veering  Patterns  In  Planetary  Gears 


When  applied  to  planetary  (or  any  epicyclic)  gears,  the  coupling  factors  reduce  to 
particularly  simple  forms  because  of  the  unique  structure  of  the  vibration  modes.  In  the 
frequency  plots  of  tuned  planetary  gears,  loci  of  the  same  type  (rotational,  translation,  and 
planet  modes)  never  cross  each  other  but  veer  away  when  they  come  close.  The  modes  of 
different  types  can  cross  each  other  but  never  switch  modes  through  veering.  This  special 
veering  pattern  is  analytically  proved  using  the  well-defined  vibration  mode  properties 
(Lin  and  Parker,  2000b).  The  general  veering/crossing  patterns  are  summarized  in  Table 
4.1.  This  pattern  is  generally  valid  when  stiffness  or  inertia  parameters  are  varied. 

The  special  veering  patterns  are  helpful  to  trace  the  evolution  of  eigenvalue  loci  and 
identify  the  effects  of  design  parameters  on  planetary  gear  vibration.  The  planetary  gear 
in  a  helicopter  powertrain  is  used  as  an  example.  The  nominal  model  parameters  are 
given  in  Table  C.l,  case  H.  Table  3.1  identifies  the  mode  type  and  where  the  dominant 
strain  energy  is  in  each  mode.  The  natural  frequencies  are  numbered  at  the  nominal 
conditions  (indicated  by  the  dashed  lines  in  Figure  4.6-1 1). 
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Table  4.1  Veering  (V)  and  crossing  (X)  patterns  in  planetary  gears.  X/V  means  crossing 
only  occurs  when  ksp=  krp  and  0(r=cy. 
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Figure  4.7  Natural  frequencies  versus  the  planet  bearing  stiffness  kp.  The  planet  mode  11 
(pure  tangential  type  PI)  crosses  the  planet  mode  15  (no  tangential  motion  type  P2).  The 
nominal  conditions  (dashed  line)  are  listed  in  Table  C.l,  case  II. 

Mesh  stiffnesses  ksp,  krp  (Figure  4.6and  4.8)  have  little  influence  on  the  low  natural 
frequencies  C0i~0)io.  This  is  because  these  modes  are  governed  by  bearing  stiffnesses 
(Table  3.1)  that  are  much  smaller  than  the  mesh  stiffnesses.  Modes  15-18  have  large 
strain  energy  in  the  sun-planet  meshes  and  are  affected  by  ksp  (Figure  4.6);  modes  1 1-14 
have  substantial  strain  energy  in  the  ring-planet  meshes  and  are  affected  by  k^,  (Figure 
4.8).  When  ksp  is  reduced  from  the  nominal  value,  the  changing  0)i5~0)i8  approach  the 
0)ji~(0i4  loci.  Because  loci  of  the  same  type  can  not  intersect,  veering  occurs  between 
rotational  modes  18  and  14,  translational  mode  pairs  (16,17)  and  (12,13),  and  planet 
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modes  15  and  11.  Below  the  veering  zones  ( ksp<100  N/jum),  modes  11-14  are  very 
similar  to  modes  15-18  above  the  veering  zones  ( ksp>800  N/jum).  In  the  same  way,  one 
can  predict  the  trend  of  frequency  loci  as  k,v  is  increased  (Figure  4.8).  Using  the  derived 
veering  patterns  and  modal  properties,  the  actual  modes  affected  by  varying  parameters 
can  be  detected  easily  although  the  plots  are  complicated  by  veering  phenomena. 

Support  stiffnesses  kh,  khll,  h=c,r,s  of  the  carrier,  ring,  and  sun  can  vary  over  a  wide 
range  depending  on  the  configuration  (fixing  or  floating  these  components).  Rotational 
and  planet  modes  are  independent  of  the  transverse  support  stiffness  kh  because  they  have 
no  translation  of  the  carrier,  ring,  and  sun;  only  translational  modes  are  affected  by 
changes  in  k/,  (Figure  4.9a, b).  Considering  the  veering  effects,  kh  significantly  affects  only 
one  pair  of  translational  modes  with  dominant  strain  energy  in  the  transverse  supports. 
When  the  rotational  support  stiffnesses  kh„  are  altered,  similar  results  are  obtained  (Figure 
4.9c, d),  except  it  is  the  rotational  modes  that  are  susceptible  to  kh„  variations. 

Altering  planet  parameters  affects  most  natural  frequencies  as  all  modes  involve 
planet  deflections,  in  general.  Applying  the  derived  veering  results,  five  pairs  of  veering 
are  identified  for  changing  planet  bearing  stiffness  in  Figure  4.7:  rotational  modes  8  and 
18,  4  and  14,  translational  modes  (9,10)  and  (16,17),  (6,7)  and  (12,13),  and  planet  modes 
5  and  15.  For  large  stiffness  kp>5000  N/fmt,  eight  natural  frequencies  increase  rapidly  to 
outside  the  range  of  interest.  Planet  mass  mp  and  moment  of  inertia  Ip  also  have 
significant  influence  on  the  natural  frequencies  (Figure  4.10). 
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4.6  Discussion  and  Summary 


Eigenvalue  loci  veering  also  occurs  when  degenerate  modes  of  symmetric  system 
are  separated  by  small  disorders.  In  planetary  gears,  the  cyclic  symmetry  can  be  broken 
by  differing  mesh  stiffnesses  at  each  planet  mesh,  manufacturing  variations,  and 
assembly  errors.  For  cyclically  symmetric  or  periodic  systems  with  small  disorders  and 
weak  structural  coupling,  mode  localization  often  accompanies  eigenvalue  loci  veering 
(Pierre,  1988).  Planetary  gears  have  relatively  strong  coupling  between  the  planets 
through  the  carrier  and  teeth  meshes,  so  strong  mode  localization  is  unlikely  even  in  the 
presence  of  loci  veering.  After  examination  of  many  cases  with  various  configurations 
and  parameters,  the  authors  have  not  found  a  realistic  example  of  mode  localization  in 
planetary  gears. 

The  special  veering  patterns  of  planetary  gear  eigenvalue  loci  are  easily 
summarized.  Two  approaching  eigenvalue  loci  of  the  same  type  (rotational  mode, 
translational  mode,  and  planet  mode)  veer  away  while  two  loci  of  different  mode  types 
cross  each  other.  The  mode  shapes  are  exchanged  across  the  veering  zone.  In  the  veering 
zone,  the  modes  are  strongly  coupled  and  markedly  different  than  outside  the  zone.  One 
can  expect  significant  differences  in  response  from  these  changed  modes.  These  rules 
result  from  planetary  gears'  unique  modal  properties  and  apply  to  all  design  parameters. 

The  effects  of  key  design  parameters  are  summarized  below: 
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1.  Mesh  stiffness.  ksp  and  kn>  each  control  three  different  natural  frequencies  associated 
with  one  rotational  mode,  one  pair  of  translational  modes,  and  one  group  of  planet 
modes.  Dominant  strain  energy  occurs  in  the  tooth  meshes  of  these  vibration  modes. 

2.  Carrier,  ring,  and  sun  parameters,  kh  and  kht„  h=c,r,s  each  affect  only  one  natural 
frequency.  The  transverse  stiffness  kh  controls  one  pairs  of  translational  modes  and  the 
torsional  stiffness  k\m  controls  one  rotational  mode.  Floating  or  fixing  the  carrier,  ring,  or 
sun  has  limited  influence  on  planetary  gear  modal  properties.  The  carrier,  ring,  and  sun 
masses  and  moments  of  inertia  affect  the  same  frequencies  as  their  corresponding  support 
stiffness,  though  the  frequencies  vary  in  the  opposite  direction. 


Figure  4.8  Natural  frequencies  versus  the  ring-planet  mesh  stiffness  krp.  Loci  12-14  will 
veer  with  loci  16-18  when  krp  is  further  increased.  The  nominal  conditions  (dashed  line) 
are  listed  in  Table  C.l,  case  n. 
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3.  Planet  parameters.  Planet  bearing  stiffness  and  planet  inertia  are  the  most  influential 
parameters  and  affect  most  natural  frequencies.  A  stiff  planet  bearing  can  be  beneficial 
for  resonance  timing  because  it  substantially  reduces  the  number  of  natural  frequencies 
in  the  lower  frequency  range  that  is  commonly  of  most  interest. 


Figure  4.9  (a),  (b)  kh,h=s,c  only  affect  translational  modes,  (c),  (d)  khu,  h=s,c  only  affect 
rotational  modes.  The  nominal  conditions  (dashed  line)  are  listed  in  Table  C.l. 
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CHAPTER  5 


PARAMETRIC  INSTABILITY  FROM  MESH  STIFFNESS  VARIATION 


In  previous  chapters,  the  mesh  stiffnesses  are  considered  constant.  Current  work 
investigates  the  impacts  of  mesh  stiffness  variation  on  planetary  gear  dynamics.  The 
variation  of  mesh  stiffness  gives  rise  to  a  parametric  excitation.  The  linear,  time-varying 
system  studied  herein  is  governed  by  a  set  of  coupled  Hill’s  equation 

Mq  +  [K0+Kv(O]q  =  F(O  (5.1) 

where  Ko  is  the  mean  stiffness  matrix  and  K v(t)  is  the  variational  part  of  the  stiffness 
matrix.  F(t)  is  the  external  excitation.  K v(t)  is  periodic  with  frequency  Q  (mesh 
frequency).  Mesh  stiffness  depends  on  element  compliance,  tooth  error,  profile 
modification,  transmitted  torque,  and  the  position  of  contact  (Kasuba  and  Evans,  1981). 
When  K v(t)  is  simple  harmonic  in  t,  (5.1)  becomes  a  set  of  coupled  Mathieu  equations. 

The  parametric  excitation  from  K v(t)  causes  instability  and  severe  vibration  when 
harmonics  of  the  excitation  frequencies  are  close  to  particular  combinations  of  the  natural 
frequencies.  Three  types  of  instability  are  of  most  interest:  (1)  primary  instability  Q~2(Op, 
(2)  secondary  instability  Q~C0p,  and  (3)  combination  instability  Q~(Op+(Oq.  The  objective 
of  this  study  is  to  systematically  determine  the  operating  conditions  leading  to  parametric 
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instabilities  in  planetary  gears.  As  discussed  in  Chapter  1,  the  existing  analysis  on 
parametric  instability  is  scarce,  inconsistent,  and  incomplete  for  multi-mesh  gears 
systems.  A  simple  two-stage  gear  chain  is  analyzed  first,  and  then  the  analysis  is 
extended  to  planetary  gears  with  multiple  meshes.  Perturbation  methods  are  used  in  this 
study  to  determine  the  boundaries  separating  the  stable  and  unstable  regions.  Floquet 
theory  and  numerical  integration  are  used  to  validate  the  analytical  findings.  The  effects 
of  contact  ratios  and  mesh  phasing  on  parametric  instability  are  quantitatively  identified. 
The  interaction  of  instability,  tooth  separation  nonlinearity,  and  dynamic  response  is 
discussed. 

5.1  Two  Stage  Gear  Systems 

5.1.1  System  Model 

Two-stage  gear  trains  have  three-gear  and  four-gear  configurations  (Figure  5.1). 
The  gear  bodies  and  intermediate  shaft  connecting  gears  2  and  4  are  assumed  rigid.  The 
tooth  meshes  are  modeled  as  linear  springs  with  stiffnesses  ku,  hi.  The  anchored  shaft  is 
flexible  with  torsional  stiffness  kio.  The  gears  have  base  radii  r„  i=l ,2,3,4.  Only 
rotational  vibrations  Oi,  62,  Os  are  considered.  The  equivalent  masses  are  mj=  Ij/rj  ,  »/?= 
h/r2  ,  and  nis=  vl '3/rf,  where  /,  are  the  moments  of  inertia  of  the  gears  and  their 
connected  shafts,  v=r/rj  for  four-gear  trains,  and  v=l  for  three-gear  chains.  The 
equivalent  stiffnesses  are  ko=  ku/rf,  kj=  ku,  and  ks=  Vku.  The  shaft/gear  rotations  are 
measured  by  the  base  radius  deflections  xj=  Ojrj,  x?=  dir 2,  and  xs=Osrs/v.  The  system 
stability  is  governed  by  the  free  vibration  equation 
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Figure  5.1  Two-stage  gear  systems  with  (a)  four  gears  and  (b)  three  gears.  kL1,  kL2 
denote  mesh  stiffnesses  and  Z2,  Z4  are  numbers  of  gear  teeth,  kio  is  the  torsional 
stiffness  of  the  anchored  shaft. 


Mq  +  [K0+Kv(O]q  =  O  (5.2) 

where  q  =  [x15  x2,x3]T.  M =diag(mi,  m2  ,  m3)  is  the  inertia  matrix.  The  stiffness  matrix  is 
represented  by  a  mean  value  Kq  and  a  variational  part  Kv(t)  as 


+ 

1 _ 

k,l 

0 

K 1 

Ki 

0 

K„  = 

kBl 

kgl+kg2 

kS2 

Kv(0  = 

K 1 

kv]  +  kv2 

K  2 

(5.3) 

0 

kg2 

kg2_ 

_o 

K  2 

K  2 . 

where  kgi  and  kvi(t)  are  the  mean  and  time-varying  components  of  the  mesh  stiffnesses, 
ki(t)  =  kgi  +  kjt)  (Figure  5.2).  The  variational  parts  Ki(t)  are  periodic  at  the  mesh 

frequency  Q,  and  expressed  in  Fourier  series  as 


NAS  A/CR— 2001-210939 


61 


/  =  1,2 


(5.4) 


MO  =  2kai^(aln  sin  sQj  +  b^  cos  sQ,t), 

/=i 

where  2kai  is  the  peak-to-peak  amplitude  of  kYi  (Figure  5.2).  The  mesh  frequencies  Oj  and 
Q2  are  related  by  0j=R02,  where  R=Z2/Z4  and  Z?,  Z4  are  the  numbers  of  teeth  on  gears  2 
and  4  (Figure  5.1a).  Note  R=l,  0]=02  for  three-gear  systems  (Figure  5.1b).  Mesh 
stiffness  variation  is  obtained  through  measurement,  calculation,  or  simple  specification 
(e.g.,  sinusoidal  or  rectangular  wave).  For  spur  gears,  rectangular  waves  are  often  used  to 
approximate  the  mesh  stiffness  alternating  between  n  and  n+1  pairs  of  teeth  in  contact 
(Kahraman  and  Blankenship,  1999a).  In  this  study,  the  kYi  are  specified  as  rectangular 
waves  with  variational  amplitudes  kai,  periods  T,  =Q/2n,  contact  ratios  c„  and  phasing 
angles  pjTj  (Figure  5.2b).  Thus, 

a,(,)  = — — sinfMc,  —  2_p,)]sin(/flc<),  b?n  = — —cos [/;r(c,  -2/?,)]sin(/;ze()  (5.5) 

In  In 

in  (5.4)  for  /  =  1, 2,  •••.  Without  loss  of  generality,  one  can  specify  pi=0,  p2=h  (h  is 

called  mesh  phasing).  In  practice,  the  first  three  or  four  Fourier  terms  reasonably 

approximate  the  mesh  stiffness  variation. 

For  the  time-invariant  case,  the  eigenvalue  problem  associated  with  (5.1)  is 
K  A  =  co? M0, .  The  vibration  modes  (j),  are  normalized  as  <bTMO=I  with  O=[0y,  (jh,  ]. 
Applying  the  modal  transformation  q  =  4>u  and  using  (5.4),  equation  (5.2)  becomes 

+££x2ei(D-  sin/Qi' + E- coslQ 10 

r=l  /= 1 

+  2s2 (F?'i  )  sin lQ2t  +  G(J?  coslQ2t)]ur  =0,  n  =  1,2,3 

where  £}  =  km  Ik  . ,  i=l,  2  and  the  matrices  D,  E,  F,  G  are 
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Figure  5.2  Modeling  of  mesh  stiffnesses  kj(t)  =  kgj  +kvi(t).  Cj  are  contact  ratios,  kgj  are 
average  mesh  stiffnesses,  and  piTj  are  phasing  angles. 
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5.1.2  Conditions  of  Parametric  Instability 

Parametric  instability  depends  on  the  frequency,  amplitude,  and  shape  of  the 
parametric  excitations.  In  gear  systems,  these  factors  are  directly  associated  with  the 
operating  speed  and  gear  design  parameters  such  as  contact  ratio,  facewidth,  diametral 
pitch,  pressure  angle,  material  properties,  and  so  on.  The  corresponding  model 
parameters  are  the  stiffness  variation  amplitudes  kai,  ka2,  mesh  frequencies  T2j,  Qi, 
contact  ratios  c;,  c?,  and  mesh  phasing  h. 

The  variation  amplitudes  kai  are  assumed  small  compared  to  the  average  mesh 
stiffnesses  kgi,  so  £,=  kjkgi  «1  in  (5.6).  At  this  point,  £;=£?=£  is  specified;  the  case  of 
£]^£~2  is  discussed  later.  Using  the  method  of  multiple  scales,  three  different  mesh 
conditions  are  examined. 

1.  Three-gear  Systems:  Equal  Mesh  Stiffness  Variations 

In  three-gear  systems  (Figure  5.1b),  the  two  meshes  have  the  same  mesh 
frequencies  Qj  =C12  =  Q.  We  consider  the  case  where  the  gear  facewidths  and  material 
properties,  which  primarily  determine  mesh  stiffness  for  a  given  number  of  teeth  in 
contact,  are  such  that  the  amplitudes  of  mesh  stiffness  variation  are  the  same  at  the  two 
meshes  (£x  =  e2  =  e).  The  contact  ratios  and  mesh  phasing  are  allowed  to  differ  between 
the  two  meshes,  however.  In  practice,  the  contact  ratios  are  changed  using  center 
distance,  diametral  pitch,  pressure  angle,  tooth  addendum,  and  other  parameters.  The 
mesh  phasing  depends  on  the  layout  of  the  gears  and  the  numbers  of  teeth. 
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The  parametric  instability  when  sC2  is  close  to  0),,+coq  is  considered.  Let 
IQ  =  (0  +  coq  +  <J£ ,  where  a  is  a  detuning  parameter.  Using  multiple  scale  method  in 


(5.6),  the  boundaries  of  the  instability  regions  are  (Lin  and  Parker,  2000c) 


a  =  %,+a,.±el A”  >  A®  =—[(02  +  Ffy  +  (£">  +  G®)2] 


For  single  mode  instabilities  (p=q ),  equation  (5.8)  becomes 


Q=TL±/^[(^>+^>)!+<^>+G")a]l'!’  ,=,A- 


As  an  example,  Figure  5.3  shows  the  boundaries  for  primary  (1=1),  secondary 
(1=2),  and  combination  (p^q,  1=1 )  instabilities.  The  parameters  are  given  in  Table  5.1  and 
Cj  =c2  =1.5,  h  =  0.  Floquet  theory,  numerical  integration  to  compute  the  fundamental 
matrix,  and  the  fundamental  matrix  eigenvalues  (Nayfeh  and  Mook,  1979)  are  used  to 
determine  the  actual  instability  regions  denoted  by  *  in  the  figures.  The  first-order 
approximations  from  (5.8)  agree  well  with  the  numerical  solution,  even  when  £  is  not 
small.  In  Figure  5.3,  the  instability  region  around  £2~2cos  is  much  larger  than  that  of  the 
primary  instabilities  around  2a>i  and  2  ah.  This  is  explained  by  examining  the  vibration 

modes.  From  (5.8),  the  primary  instability  boundary  slopes  are  governed  by  ^A^". 
Expansion  of  D,  E,  F,  G  in  (5.7)  yields  the  diagonal  terms 


+A,)je  e »> = ktMP + a,)j^ 
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where  fop,  fop,  fop  are  the  gear  rotations  in  mode  fo.  Let  Sip=fop+fop  =xj+X7  be  the 
relative  deflection  of  the  first  mesh  in  mode  fo.  Similarly,  fop-fop+fop  =x?+X3  represents 
the  modal  deflection  in  the  second  mesh.  For  the  primary  instability  boundary  around 
2cOp,  insertion  of  (5.10)  into  (5.8)  yields 

a'"  =  K*,AX’ +W)1  +<*,AX’  +kJlrb (5.i  i) 

The  mesh  deflections  Si,  fo  in  each  mode  can  be  observed  from  the  mode  shapes  (Figure 
5.4).  The  two  meshes  in  fo  are  both  in  phase  and  have  smaller  Si,  fo  than  those  of  fot 
where  the  two  meshes  are  both  out  of  phase.  Thus,  A(^  >  A*1]  and  the  instability 

boundaries  around  2a ^  have  larger  slope  than  those  around  2coj.  Mode  fo  has  one  mesh  in 
phase  and  the  other  out  of  phase,  so  the  size  of  its  primary  instability  regions  is  between 
that  of  fo  and  fo  ( A33  >  A(^  >  A*,1} ).  In  addition,  the  mesh  deflections  in  a  vibration  mode 

are  related  to  the  modal  strain  energy  Ul  =  kx8x  12,  U2  =  k2S;  / 2 .  Examination  of  (5.1 1) 
shows  that  vibration  modes  with  more  strain  energy  in  the  meshes  have  larger  instability 
regions  and  are  more  susceptible  to  parametric  excitations.  The  above  results  apply  for 
mesh  stiffness  variations  of  arbitrary  shape. 

For  mesh  stiffnesses  having  rectangular  waveforms,  one  can  clearly  identify  the 
effects  of  contact  ratios  and  mesh  phasing  on  the  instability  regions.  Use  of  (5.5)  in  (5.1 1) 
yields 

Ap],  =  [k2glS ,p  sin2  {cx7t)  +  k]2S2p  sin \c2n)  + 

,  ,  2  ,  (5-12) 

2k  xk  28lp82  sin(c^)  sin(c2^r)  cos^  -c2+  2h)n]{ - ) 

na)„ 
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Inertia 

mj=l,  1)12=0.3,  ms=4.0 

Average  mesh  stiffness 

II 

Shaft  stiffness 

ko=0.5 

Contact  ratio 

1<  Cl,  c2<2 

Mesh  phasing 

pi=0,  0<p2~h<l 

Variational  amplitude 

0<  kai,  ka2  <0.5 

Table  5.1  Parameters  of  an  example  system  in  Figure  5.2. 


Figure  5.3  Instabilities  regions  when  Q.l=Q2=Q,£l=£2=£;  —  analytical  solution; 
***  numerical  solution.  The  parameters  are  from  Table  5.1  and  cx  -  c2  =  1.5,  h  =  0. 
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The  minimum  value  A®  =  0  is  obtained  when  sinc^  =  sin  c27l  =  0  and  the  unstable  region 

vanishes.  This  is  achieved  for  integer  contact  ratios  cj,  c2  where  the  number  of  tooth  pairs  in 
contact  remains  constant  and  mesh  stiffnesses  are  time-invariant.  For  given,  non-integer  c12 

between  1  and  2,  is  minimized  by  setting  cos(c,  -  c2  +  2h~)K  =  -1  or  c,  -  c2  +  2h  =  ±1, 3  . 

By  properly  choosing  the  contact  ratios  and  mesh  phasing,  the  parametric  instability  regions  can 
be  dramatically  reduced.  Figure  5.5  compares  the  instability  regions  for  three  cases.  The  most 
severe  condition  for  primary  instabilities  (dash-dot  lines)  is  c1=c2  =  1.5  and  h  =  0,  which 

maximizes  A(*j,  in  (5.12).  This  condition  is  markedly  improved  by  changing  the  phasing  h  =  0.5 
so  that  Cj  —  c2  +  2h  =  1  (dashed  lines).  When  the  contact  ratios  are  close  to  integers  (solid  lines, 
c,  =  1.1,  c2=  1.9,  h  =  0.4 ),  the  primary  instability  region  becomes  even  smaller. 

Similar  conditions  are  obtained  for  secondary  and  combination  instabilities  (Lin 
and  Parker,  2000c).  Unfortunately,  the  primary,  secondary,  and  combination  instability 
regions  cannot  be  minimized  at  the  same  time.  The  conditions  reducing  the  primary 
instability  regions  (dashed  lines,  Figure  5.5)  result  in  large  combination  instability 
regions,  and  vice  versa  (dash-dot  lines).  Depending  on  specific  applications,  a  trade-off 
may  be  made  to  reduce  multiple  instability  regions,  though  none  are  true  minima  (solid 
lines).  Adjusting  contact  ratios  and  mesh  phasing  is  clearly  an  effective  means  to 
minimize  instability  regions  and  avoid  resonances  under  operating  conditions. 

2.  Four-gear  Systems:  Equal  Mesh  Stiffness  Variations 

Two-stage  countershaft  systems  (Figure  5.1a)  have  two  different  mesh  frequencies  £2i=RQ2, 
which  means  more  instability  regions  than  three-gear  systems.  We  consider  the  case  where  the 
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gear  facewidth  and  material  are  such  that  the  mesh  stiffness  amplitudes  are  identical  at  the  two 
meshes  (eI=e2=e),  although  the  contact  ratios  and  phasing  are  not  restricted.  Depending  on  the 
ratio  R=Z2/Z4,  the  parametric  instability  regions  associated  with  12/  and  Qi  may  overlap  each 
other.  For  R  =  m/j  (m,  j  are  integers),  the  l=j  instabilities  (single  mode  and  combination)  of 
and  the  l-m  instabilities  of  12?  occur  simultaneously.  Because  their  instability  regions  are 
typically  the  largest,  the  interactions  involving  either  m  or  j=l  are  of  most  interest. 


Figure  5.4  Vibration  modes  for  the  time-invariant  system  with  parameters  in  Table  5.1. 
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Figure  5.5  Comparison  of  instability  regions  for  various  contact  ratios  and  mesh  phasing. 
The  parameters  are  in  Table  5.1.  —  cx  =c2  =1.5,  h  =  0.5; — c,  =1.1,  c2  =  1.9,  h  =  0.4;  - 
--  q  =c2  =1.5,  h  =  0. 


When  i?  ^  m,  1/m  for  integer  m,  the  1=1  instabilities  from  one  mesh  decouple  from 
the  l=m  instabilities  of  the  other  mesh.  In  this  case,  instability  occurs  when  IQj  or  IQo  is 
close  to  COp+COq ,  but  these  instability  boundaries  can  be  calculated  independently.  For  IQi 
=  Cty+Cty+eGi,  the  condition  separating  bounded  and  unbounded  solutions  is 


Q,  = 

Q2  = 


s\n(lcx7r) 


l 


<°p% 


(5.13) 
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The  stability  regions  associated  with  each  mesh  frequency  depend  on  the  individual 
contact  ratios  but  are  independent  of  mesh  phasing  as  the  two  mesh  excitations  are 
uncoupled.  For  primary  and  combination  instabilities  (1=1),  maximum  regions  occur 
when  c\  =  c2  =  1 .5  and  minimum  regions  require  Cj=c2={1,  2}.  For  secondary 
instabilities  (1=2),  the  maximum  and  minimum  conditions  are  c,  =  c2  ={1.25,  1.75}  and 
c\  =c2  ={1,  1.5,  2},  respectively.  Figure  5.6a  shows  the  instability  regions  for  R=3/5. 
The  primary  instability  associated  with  C22=2(Q3+£<?2  occurs  at  Q\  =R£22=R(2c03+e<72)  in 
Figure  5.6a. 

When  R  =  m  or  1/m  for  integer  m,  the  parametric  excitations  from  the  two  meshes 
interact.  Consider  the  case  with  R=l/m,  where  the  l=m  instabilities  caused  by  Qj  overlap 
with  the  primary  instabilities  caused  by  Qi-  Considering  instability  of  the  p- th  mode 
where  £2i=2(Op/m+eai  and  £22=mQi=2cOp+mecTi,  the  boundaries  have 

^  =  [(D™  +02  +(£«  +G^)2f2  !(mcop)  (5.14) 

For  example,  when  TT=2i2/  (R=l/2),  the  boundaries  for  Qi  secondary  instabilities 
(overlapping  with  ^primary  instabilities)  are 

a,  =<»p±^-l(D%+FVy+(E%+G%)2]'n  (5.15) 

Figure  5.6b  shows  instability  regions  in  the  (12/,  s)  plane  for  R=l/2.  Note  the  instability 
at  £2]~C0s  couples  with  the  instability  at  &2~2aj3,  and  the  combined  instability  region  is 
much  larger  than  the  case  without  interaction  (Figure  5.6a).  Using  (5.5)  and  (5.10)  in 
(5.15),  the  slopes  of  these  boundaries  are 
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(5.16) 


<JX  =  [klA*P  sin2(2c1^)  +  k2g282p  sin2(c,;r) 

+  2 kglkg2SXpS2p  sin(2c,;r)  sin(c2;r)  cos(2c,  -c2+  2h)tt]1'2  /  cop 

Minimization  of  0>  requires  c,={l,  1.5, 2}  and  c2={l,  2}  for  1<c12<2.  The 

instability  region  can  also  be  reduced  by  adjusting  the  phasing  h  according  to 
cos(2cj  -  c2  +  2 h)n  =  ±1  with  sign  the  same  as  sin(2c,;r) .  The  primary  instability  regions 
under  Qj  (1=1)  do  not  coincide  with  any  other  instability  regions.  Other  overlap  situations 
are  possible,  such  as  the  22?  secondary  instability  (1=2)  overlaps  with  the  Qi  fourth 
instability  (1=4),  but  the  interaction  between  these  higher  instabilities  is  typically  weak 
and  the  instablity  regions  are  much  smaller.  Combination  instabilities  can  be  analyzed 
similarly. 

3.  Three  and  Four-gear  Systems:  Unequal  Mesh  Stiffness  Variations 

This  general  case  allows  all  parameters  of  the  two  mesh  stiffnesses  to  differ.  In 
contrast  with  prior  cases,  the  gears  may  have  differing  facewidths  and  material  properties 
such  that  the  amplitudes  of  stiffness  variation  at  each  mesh  vary  independently  (£> **£?)• 
The  contact  ratios  and  mesh  phasing  are  unrestricted.  The  design  of  one  mesh  must 
account  for  dynamic  interactions  with  the  mesh  stiffness  variation  of  the  other. 

When  R  #  m,  1/m  for  integer  m,  there  is  no  interaction  between  the  parametric  excitations 
from  the  two  meshes.  The  Qi  instabilities  are  only  affected  by  £j  and  the  instabilities 
are  only  affected  by  £?. 

When  R  =  m  or  1/m  for  integer  m,  a  mode  may  be  simultaneously  driven  to 
instability  by  both  mesh  excitations.  In  this  case,  the  first  mesh  instability  regions  can  be 
significantly  affected  by  the  presence  of  £?  and  vice  versa.  Closed-form  boundaries  of  the 
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form  (for  primary  instability)  Oi  =  2cq}+£i<j]+£2<7?  for  independently  varying  £j,  £2  are 
cumbersome.  Alternatively,  simple  yet  accurate  approximations  for  the  instability 
boundaries  are  obtained  by  presuming  a  linear  variation  of  the  boundaries  in  the  £1) 
plane  for  given  £?.  To  construct  this  linear  approximation,  one  point  is  calculated  under 
the  condition  £l=0,£2=C  and  a  second  point  is  obtained  at  £l=£2=C .  The  primary 
stability  boundary  limits  for  £x  =0,£2  =  C  are 

£2,  =  2®,  ±— [(F«)2  +(G™)!],,!  (5.17) 

From  (5.9),  the  stability  boundary  limits  for  £x  =  C,e2  =  C  are 

fl,  =2®,  ±— [(Z>«  +  F«>)2  +(£“  +G®)2]1'2  (5.18) 

An  example  is  for  the  primary  instability  when  R =7,  C=0.3.  Connecting  the  two 
points  obtained  from  (5.17)  and  (5.18)  yields  the  instability  boundaries,  which  agree  well 
with  the  numerical  solution  (Figure  5.7a).  Assembling  the  (£2i,  £})  planes  for  various  £o=C 
generates  three-dimensional  plots  of  versus  £1,  £2  (Figure  5.7b).  The  parametric 
excitation  in  the  second  mesh  dramatically  changes  the  shapes  of  the  instability  regions. 
Notice  that  the  second  parametric  excitation  widens  the  primary  instability  region  for 
small  £1  compared  to  mono-frequency  excitation  (Figure  5.3).  In  contrast,  the 
combination  instability  at  Qi~Ob  +a>3  disappears  near  £j=0.23  in  Figure  5.7a.  In  other 
words,  an  otherwise  unstable  system  is  stabilized  by  the  presence  of  a  second  parametric 
excitation.  The  solid  line  in  Figure  5.7b  indicates  points  where  the  Ob  +CO3  combination 
instability  vanishes. 


NAS  A/CR— 2001-210939 


73 


(a) 


3.5 

3 

2.5 1 - 1 - 1 - 1 - 1 - 1 - 1 

0  0.05  0.1  0.15  0.2  0.25  0.3 


2Ro)„ 


=T 


j 


€ 


Figure  5.6  Instabilities  regions  when  Q,  =  RQ,2,  £:  =£-,=£ .  (a)  R=3/5,  (b)  R=l/2.  The 
parameters  are  in  Table  5.1  and  c,  =  c2  =  1.5,  h  =  0 .  ***  denotes  numerical  solutions. 
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and  the  solid  line  indicates  vanishing  of  the  combination  instability.  The  parameters  are 


in  Table  5.1  and  cl=c2=  1.5,  h  =  0 . 
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5.1.3  An  Example 

The  two-stage  gear  system  (Figure  5.1a)  studied  by  Tordion  and  Gauvin  (1977)  and 
Benton  and  Seireg  (1980b)  is  used  as  an  example.  These  two  papers  come  to  markedly 
different  conclusions  as  discussed  below.  The  system  parameters  are  given  in  Table  5.1 
and  Cj  =1.47,  c2  =1.57 .  In  keeping  with  the  published  work,  the  double-tooth  contact 
mesh  stiffness  kmax=l  is  kept  constant,  and  the  average  mesh  stiffnesses  kgi  decrease  as  kai 
is  increased  (Figure  5.2b). 

Tordion  and  Gauvin  assumed  that  kvi  and  kv?  have  the  same  amplitude  and 
frequency  but  different  contact  ratios  and  phasing.  They  applied  an  infinite  determinant 
method  (Bolotin,  1964)  to  plot  the  boundaries  of  primary  and  secondary  instabilities 
(dashed  lines  in  Figure  5.8).  Their  results  deviate  significantly  from  the  numerical 
solution  as  a  result  of  analytical  errors.  In  addition,  the  Fourier  expansion  they  derived  for 
rectangular  waveforms  (equations  (11)  and  (12)  in  Tordion  and  Gauvin,  1977)  is 
incorrect.  Nevertheless,  they  conclude  that  “The  phase  displacement  between  the 
meshing  stiffnesses  has  a  great  influence  on  the  width  of  the  instability  regions.” 

Benton  and  Seireg  (1980b)  considered  the  same  system.  They  decoupled  the 
equations  using  the  modal  transformation  and  neglected  the  off-diagonal  terms  of  the 
transformed  time-varying  stiffness  matrix  (that  is,  OrKv(t)<I>  )•  These  treatments  reduce 

(5.6)  to  three  uncoupled  Mathieu  equations.  The  average  value  of  two  contact  ratios  was 
used  to  make  the  stiffness  variations  kvj  and  identical.  With  these  approximations, 
they  conclude  that  the  instability  regions  are  independent  of  the  mesh  phasing,  that  is, 
“the  normal  mode  technique  . . .  without  considering  the  phase  variations  . . .  provide(s)  a 
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relatively  simple  means  of  predicting  the  instability  regions  with  sufficient  accuracy  for 
practical  purposes.”  This  conflicts  directly  with  Tordion  and  Gauvin.  In  fact,  the  mode 
uncoupling  method  does  not  provide  satisfactory  results  when  the  mesh  phasing  is  non¬ 
zero  (Figure  5.8b). 

The  perturbation  results  resolve  the  discrepancy:  Mesh  phasing  strongly  impacts  the 
mesh  stiffness  variation  instabilities.  The  excellent  agreement  of  analytical  and  numerical 
stability  boundaries  confirms  this  finding  (Figures  5. 3-5. 8). 

To  further  validate  the  stability  conditions,  free  responses  under  non-trivial  initial 
conditions  are  calculated  numerically  (Figure  5.9)  for  the  parameters  at  point  A  of  Figure 
5.8  (12  =  4.2,  ka=e=0.3).  For  point  A  in  Figure  5.8a,  the  responses  are  unstable  (Figure 
5.9a),  as  identified  by  perturbation  and  numerical  methods.  This  point,  however,  is  stable 
according  to  Tordion  and  Gauvin  (Figure  5.8a).  When  the  phasing  h=0.4  at  point  A 
(Figure  5.8b),  stable  responses  occur  (Figure  5.9b).  This  is  consistent  with  the 
perturbation  and  numerical  solutions  but  conflicts  with  both  Tordion  and  Gauvin’s  and 
Benton  and  Seireg’s  results. 

5.1.4  Discussion 

Rectangular  waveforms  are  close  approximations  of  the  mesh  stiffness  in  spur  gears 
with  involute  teeth.  For  helical  gears  or  spur  gears  with  tooth  modification,  mesh  stiffness 
deviates  from  the  rectangular  shape.  Equation  (5.5)  is  not  valid  for  other  functions,  but 
the  general  Fourier  expansion  (5.6)  can  still  be  used  in  matrices  D,  E,  F,  G  to  determine 
the  instability  boundaries. 
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If  damping  is  considered,  the  system  stability  improves  and  the  instability  regions 
shift  to  the  right  in  the  (i2,  e)  plane.  Furthermore,  damping  and  non-linearity  (e.g.,  tooth 
separation)  must  be  considered  to  determine  the  limit  cycle  amplitude  of  the  dynamic 
response  when  operating  conditions  cause  instability. 

The  instability  analysis  can  be  reduced  to  single  mesh  gears  with  one  natural 
frequency  co,,.  From  (5.9),  primary  and  secondary  instabilities  vanish  as  the  contact  ratio 
cj={1,2}  and  C]={1,  1.5,  2},  respectively.  Maximum  primary  and  secondary  instability 
occurs  at  C]=1.5  and  ci~{1.25,  1.75},  respectively.  Kahraman  and  Blankenship  (1999) 
experimentally  studied  a  pair  of  spur  gears  under  mesh  stiffness  excitation  for  various 
contact  ratios.  They  showed  that  the  amplitude  A  /  in  the  first  mesh  frequency  harmonic  of 
the  response  is  minimized  when  the  contact  ratio  C]={1.0,  2.0}.  This  is  because 
parametric  excitations  are  eliminated  for  integer  contact  ratios.  When  the  mesh  frequency 
C2~(On,  their  measured  A /  reaches  maximum  for  ci~1.4.  A  possible  explanation  for  the 
high  A i  at  this  contact  ratio  is  due  to  the  combined  effects  of  primary  and  secondary 
instabilities.  First,  for  both  primary  instability  excited  by  the  first  harmonic  of  k(t)  and 
secondary  instability  excited  by  the  second  harmonic  of  k(t),  the  dominant  response 
frequency  is  (On.  We  now  examine  the  contact  ratios  where  both  instabilities  are  active. 
Although  the  maximum  primary  instability  region  occurs  at  C]=1.5,  the  secondary 
instability  region  is  eliminated  there.  For  ci=1.25,  the  secondary  instability  region  is 
maximal  but  the  primary  instability  region  is  small.  For  cj~1.4  (average  of  1.25  and  1.5) 
or  1.6  (average  of  1.75  and  1.5),  however,  both  primary  and  secondary  instabilities  have 
significant  instability  regions.  Generally,  the  larger  an  instability  region,  the  higher 
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response  amplitude  occurs  due  to  this  instability  (to  see  this  heuristically,  note  that  both 
the  slope  of  the  stability  boundaries  in  (17)  and  excitation  of  first  order  response  in  (10) 
are  proportional  to  the  same  quantities  D,  E,  F,  and  G).  Accordingly,  for  ci~1.4  and  1.6, 
both  instabilities  induce  large  response  and  jointly  contribute  to  large  A /.  When  the  mesh 
frequency  i Q~2(On  (not  shown  in  Kahraman  and  Blankenship,  1999),  the  instability  is 
caused  only  by  primary  instability  excited  by  the  first  harmonic  of  k(t).  Because  the 
primary  instability  region  is  maximal  at  cj =1.5,  the  response  amplitude  A  j  around  Q~2co„ 
also  becomes  maximal.  Therefore,  from  the  viewpoint  of  dynamic  instability  and 
amplitude,  contact  ratios  in  the  range  1.4-1. 6  are  harmful  to  single-mesh  gears  at  high 
speeds. 

5.2  Planetary  Gear  Parametric  Instability 

Planetary  gears  have  multiple  time-varying  mesh  stiffnesses  and  their  parametric 
instability  has  not  been  analytically  investigated.  In  planetary  gears,  parametric 
excitations  are  complicated  as  different  contact  ratios  and  phasing  conditions  exist 
between  the  sun-planet  and  the  ring-planet  meshes.  The  objective  of  this  study  is  to 
analyze  parametric  instability  excited  by  different  time-varying  mesh  stiffnesses  in 
planetary  gears.  The  perturbation  method  used  for  two-stage  gears  is  extended  to 
planetary  gears  that  have  degenerate  vibration  modes  due  to  their  cyclic  symmetry.  The 
well-defined  modal  properties  of  planetary  gears  are  used  to  derive  simple  expressions 
for  instability  boundaries  separating  the  stable  and  unstable  regions.  From  these 
expressions,  the  effects  of  contact  ratios  and  mesh  phasing  are  analytically  determined  to 
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Figure  5.8  Comparison  of  instability  regions.  The  parameters  are  from  Table  5.1, 
Cj=1.47,  c2=1.57,  and  phasing  (a)  h=0,  (b)  h=0.4.  —  Perturbation  method;  *** 
Numerical  method;  —  Tordion  and  Gauvin  (1977);  —  Benton  and  Seireg  (1980b). 
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provide  insight  into  planetary  gear  design.  In  practice,  planet  mesh-phasing  schemes  are 
often  applied  to  cancel  or  neutralize  resonant  response  at  speeds  where  the  mesh 
frequency  is  near  a  natural  frequency  (Seager,  1975;  Kahraman,  1994b;  Kahraman  and 
Blankenship,  1994;  Parker,  2000).  In  this  same  spirit,  this  study  shows  that  particular 
parametric  instabilities  can  be  eliminated  under  certain  phasing  conditions  that  can  be 
achieved  by  proper  selection  of  design  parameters.  Tooth  separation  nonlinearity  induced 
by  parametric  instability  is  numerically  simulated. 


5.2.1  System  Model  and  Mode  Properties 

The  planetary  gear  dynamic  model  used  is  based  on  the  one  developed  in  Chapter  2 
(Figure  2.2).  Translational  degrees  of  freedom  in  that  model  are  eliminated,  and  only 
rotational  motions  of  the  gear  bodies  are  considered.  The  sun-planet  and  ring-planet  tooth 
meshes  are  modeled  as  linear  springs  with  time-varying  stiffnesses  ksn(t),  kru(t),  n=l,...N. 
The  system  equations  of  motion  are 


Mq  +  K(/)q  =  F 

M  =  diag[lc  !r]  mp,  I r  / r2,  Is  / r2 ,  I  p  / rp,  ■  ■  • ,  I p  / r2p  ] 


-IX, 

K\-K i  •" 

krN-k; 

I  *„ 

0 

~kn  "• 

~  k,N 

••• 

k,N 

K(/)  = 

K\+ks i  0 

...  0 

0 

symmetric  krN  +  ksN 


(5.19) 


(5.20) 
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q  =  [uc,ur,us,uv---,uNJ ,  F  =  [Te/re,  Tr/rr,  TJrs,  0,— ,0f 


where  the  summation  index  n  ranges  from  1  to  N.  In  the  stiffness  matrix  K (t),  each  mesh 
stiffiiess  is  represented  by  kjt)  =  ksp  +  kjt) ,  kjt)  =  krp+  k2n(t),  n-1,  ...N,  where  ksp, 

krp  are  mean  values  and  ki„,  k2„  are  time-varying  components  of  the  n- th  sun-planet  and 
ring-planet  meshes.  For  spur  gears,  rectangular  waves  are  often  used  to  approximate 
mesh  stiffnesses  alternating  between  d  and  d+1  pairs  of  teeth  in  contact.  Figure  5.10 
shows  the  mesh  stiffness  variations  ki„,  k^„  with  peak-to-peak  amplitudes  2kvs,  2kvr, 
contact  ratios  cs,  cr,  and  phasing  angles  y„T,  (ySr+Y™)T.  The  sun-planet  and  ring-planet 
meshes  have  the  identical  mesh  frequency  Q=2tt2T.  The  sun-planet  mesh  stiffnesses 
between  planets  differ  only  by  a  time  transition  (or  phase  angle),  and  likewise  for  the 
ring-planet  meshes.  Mesh  phasing  in  planetary  gears  is  determined  by  planet  position 
angles  y/„  and  the  numbers  of  teeth  zs,  zr,  zp  for  the  sun,  ring,  and  planets  (Kahraman  and 
Blankenship,  1994).  ysn=\f/^zs  /(2k)  denotes  the  mesh  phasing  between  the  first  and  n- th 
sun-planet  meshes  (ysi=0);  yrn=\f/^zr/(2K)  is  the  mesh  phasing  between  the  first  and  n-th 
ring-planet  meshes  ( yri=0 );  ysr=zp/2  is  the  mesh  phasing  between  the  sun-planet  and  ring- 
planet  meshes  for  each  planet.  Note  ysr=0  for  even  zp  and  ysr=l/2  for  odd  zp.  Expansion  of 
kj„,  k2„  in  Fourier  series  yields 


K  (0  =  2*„  X  (a®  sin  Ku  +  cos  ICU), 

!=\ 

k2„  (0  =  2kw  J  (a™  sin  l£lt  +  cos  lQ.t) 

1=1 


(5.21) 


1 
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2  2 

=-— sint/^c,  -2^„)]sin(/^),  =  — —  sin[/^(c,.  -2yir  -2/,7i )]sin(/;zu,.) 

//r  In 

2  2 

b(s'J  =  ——cos[ln(cs-2ysn)]sm(lncs),  b%  =  cos[/^(c,.  -2^,.  -2ym)]sin(/mrr) 

/;r  In 

Mesh  stiffnesses  depend  on  many  parameters  including  the  number  of  teeth  in  contact, 

gear  facewidth,  material  properties,  profile  modifications,  and  applied  load  (Kasuba  and 

Evans,  1981).  Let  £y  =krs/k  and  e2  =  kvr/k  be  the  relative  amplitudes  of  mesh 

stiffness  variation.  In  the  simplest  approximation,  mesh  stiffnesses  are  assumed 
proportional  to  the  number  of  tooth  pairs  in  contact,  that  is,  ksp=cs  kjh  krp=c,  k3,  and  kvs~ 
kit/2,  kw=  k3l/2,  where  kn,  k3,  are  one-pair  tooth  bending  stiffnesses  of  the  sun-planet  and 
ring-planet  meshes.  With  this  simplifying  assumption,  the  nominal  amplitudes  of  stiffness 
variations  are  £i=1/(2cs),  e3=l /(2cr),  and  one  obtains  the  explicit  relation  £3=£i  cs/cr.  In 
practice  and  in  the  analysis  that  follows,  however,  £j  and  e3  are  not  constrained  to  these 
“nominal”  values  but  vary  independently  of  the  contact  ratios  to  account  for  the  many 
factors  influencing  mesh  stiffness  variation  amplitudes.  For  this  modeling,  we  let 
£i=z£o/g=£  where  g=0(l)  (and  g=cs/cr  under  the  simplifying  assumption  noted  above). 
Substitution  of  (5.4)  and  (5.21)  into  (5.20)  yields 

K(0  =  K0  +  2eJ](Kvi)  sin  IQt  +  K$  cos/Qt)  (5.22) 

/=i 

where  the  average  stiffness  matrix  Kohas  the  same  form  as  (5.20)  with  ks„,  kn,  substituted 
by  ksp,  krp.  The  Fourier  coefficient  matrices  are  also  in  the  form  (5.20)  with  ks„,  kn, 
substituted  by  kvai2,gkn>a<£  for  K^and  substituted  by  k^^gk^  for  K^. 

For  the  time-invariant  case,  the  eigenvalue  problem  associated  with  (5.19)  is 
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Ko0,-  =  0)?M.<f>( 


(5.23) 


The  natural  frequencies  and  vibration  modes  have  unique  properties  (3.7),  (3.14),  (3.20). 

We  specialize  these  properties  for  the  case  of  a  rotational  vibration  model  with  fixed  ring 

and  L=N+2  degrees  of  freedom.  In  this  case,  all  vibration  modes  can  be  classified  into 

one  of  three  categories:  (1)  a  rigid  body  mode  ( a>i=0 ),  (2)  two  modes  with  distinct  natural 

frequencies  (0)2,  (0[),  and  (3)  a  group  of  degenerate  modes  with  multiplicity  N-l 

( a>3= ...  =(On+i).  In  the  distinct  modes,  all  planets  have  identical  motion 

u„=  uj,  n=l,2,...,N  (5.24) 

In  the  degenerate  modes,  the  carrier,  ring,  and  sun  have  no  motion,  and  the  planet 

rotations  satisfy 

un  =[u1sm(y/2-y/n)  +  u2siny/n]/sm\ff2,  n=l,2,...,N  (5.25) 

These  well-defined  properties  are  valid  not  only  for  equally  spaced  planets  with  position 
angles  y/„  =2n(n-l)/N,  but  also  for  diametrically  opposed  planets  with  %+n/2  =  %  +X- 

5.2.2General  Expression  for  Instability  Boundaries 

We  determine  the  operating  conditions  (that  is,  mesh  frequency  Q  and  stiffness  variation 
amplitude  e)  that  lead  to  instability  when  IQ-cOp+co^  for  integer  /.  The  rigid  body  mode 
(0)1=0)  is  not  excited  under  operating  conditions,  and  does  not  affect  the  instabilities  of 
other  modes.  Only  the  two  distinct  modes  and  the  group  of  degenerate  modes  are 
considered  in  what  follows.  The  parametric  instability  when  IQ  is  close  to  QJp+CDq  is 
considered.  Let  IQ.  =  cop  +  coq  +  ea ,  where  a  is  a  detuning  parameter  to  be  determined. 

From  (5.8),  when  (ty  and  C0q  are  both  distinct,  the  instability  boundaries  are 
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Figure  5.10  Modeling  of  mesh  stiffness  variations  in  the  sun-planet  and  ring-planet 
meshes.  cs,  cr  are  contact  ratios,  and  yn, ,  ysr  are  mesh  phasing. 


single  mode  instability:  Q  =  ± 

(5.26) 

coD  +  (Oa  £  j  777 - 

combination  instability:  Q  =  — — — —  ±  —  /{copcoq ) 

(5.27) 

where  =(£>“)’+ (O’ • 

When  (Op  or  (Oq  are  degenerate,  the  complexity  of  instability  solutions  depends  on 
the  multiplicity  of  the  degenerate  natural  frequencies.  We  first  study  the  case  with 
multiplicity  two,  say  ty,=(Or  The  single-mode  instability  boundaries  for  6^=6%  are  (Lin 
and  Parker,  2000d) 
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(5.28) 


Q=^±.f 


l  lcor 

p 

rw  =  +£<%’ +|[(<’):  +(fl"’)2 +(£"’)’  +U©9 


When  (0j,=0)q  are  degenerate  but  is  distinct,  the  combination  instability  is 


Q  =  £t±£i±£.  +  DmDm  +e^e(^  +  E^E^)/{mrcor)}vl  (5.29) 


When  the  degenerate  natural  frequencies  C03=...=0)m+2  have  multiplicity  m  >  2,  the 
critical  <ris  obtained  by  requiring  the  real  parts  of  the  eigenvalues  of  the  following  matrix 
to  be  non-positive, 

+  (<JU>3  )I„IX„,  -E«  J 

where  D^,E^  are  mxm  submatrices  of  D^,  E®  including  rows  and  columns  3,  .... 

m+2.  Generally,  no  closed-form  solution  can  be  derived  for  these  single-mode 
instabilities,  but  numerical  evaluation  can  determine  <7  and  the  instability  boundaries 
unless  D^j,E^  have  special  features  (e.g.,  diagonal  matrices)  as  discussed  later.  More 
is  possible  for  combination  instabilities.  For  combination  instability  of  distinct  (Or  and 
degenerate  (03=... =(Om+2,  the  instability  boundaries  are 

+  EZ>E%-)Ha>,ah)  (5-31) 

*  *  q=3 

The  above  expressions  apply  for  a  general  system  with  degenerate  natural 
frequencies.  These  results  reduce  to  simple  forms  when  specialized  to  planetary  gears. 
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5.2.3.  Planetary  Gear  Parametric  Instability 

Planetary  gears’  well-defined  modal  properties  are  now  used  to  simplify  the  above 
instability  conditions  to  compact  expressions  suitable  for  use  in  applications.  Mesh 
stiffness  are  approximated  by  rectangular  waveforms  (5.21).  Expansion  of  D^1,  Effi  gives 


«=1  M  =  1 

£(,)=0rK(,?V  =k  YbV)8{p)8{q)  +  k  gfb^S^S^ 

pq  rp  v2rq  sd  /  jsti  sti  sn  rpo  ni  m  m 


(5.32) 


where  8[p)  =  us  -  uc  +  un  is  the  deformation  of  the  w-th  sun-planet  mesh  in  mode  (j)p  and 


SY  =  —uc  -  un  is  the  defonuation  of  the  /7-th  ring-planet  mesh.  For  different  spacing  and 

phasing  of  planets,  matrices  D^,  E®  have  special  features  that  simplify  the  instability 
conditions. 


1.  Equally  spaced  planets 

For  equally  spaced  planets  (y/„  =2n(n-l)/N),  the  vibration  modes  have  structured 
properties  (5.24)  and  (5.25).  In  this  case,  (zs+zr)/N  =  integer  and  the  planet  meshes  have 
phasing  ys„=(n-l)zs/N  and  ym=(n-l)zr/N,  n=l,...,N.  The  two  possible  phasing  conditions 
are  examined  below. 


(1)  In-phase  Planet  Meshes 

Consider  the  case  when  all  the  sun-planet  meshes  are  in-phase  {ys„=0,  n=l,...,N) 
and  all  the  ring-planet  meshes  are  in-phase  (ym=0,  n=l,...,N).  This  design  is  typical  when 
optimizing  load  sharing  and  results  when  zs  and  zr  are  each  integer  multiples  of  N. 
However,  there  is  a  constant  phasing  ysr=0  (even  zp)  or  1/2  (odd  zp)  between  the  sun- 


NASA/CR— 2001-210939 


88 


planet  and  ring-planet  meshes  for  each  planet.  The  Fourier  coefficients  a(s'J ,a^  ,b(JJ ,b(^ 
in  (5.21)  are  independent  of  the  planet  index  n. 

For  a  distinct  natural  frequency  tq,,  vibration  mode  property  (5.24)  leads  to 
S(J]  =  =  S(r f  for  any  planet  n.  Using  (5.21)  and  (5.32)  on  (5.26),  the  primary 

instability  boundaries  have 

A®  =  (— )2[(^rt)4(*¥sin<r,)2  +  (.Sf)\k,rgsmmrf  + 

7t  (5.33) 

2  Kp^giS^S^)2  sm(7tcs )  sin(^cr )  cos  n(cs  -cr  +  2ysr)\ 

Obviously,  A(^,=  0  when  cs,  cr  are  integers  and  all  instabilities  vanish.  If  cs,  cr  ^  integer,  a 

second  choice  to  reduce  the  instability  regions  is  to  set  the  third  term  in  (5.33)  to  be 
negative  by  adjusting  cs,  cr,  and  ysr.  From  (5.27),  the  combination  instability  boundaries 
for  two  distinct  modes  have 

A(Ji  =  sin ncsf  +{k  sin ncrf  + 

K  (5.34) 

2  KpkrpS8(fSliSn)Slf  sin(^“s.)sin(^,.)cos^(cs  -cr  +  2  yj] 


Inertias  (kg) 

Ic/rea=6,  Is/rs2=2.5,  lp/rp2=2 

Planet  mass  (kg) 

mp=4 

Mesh  stiffness  (N/m) 

°°o 

II 

II 

g. 

Natural  frequencies 
(kHz) 

CQl=0,  0)2=1.212, 

0)3=ca/=1.592,  0)5=2.196 

Table  5.2  Parameters  and  natural  frequencies  of  an  example  planetary  gear  with  fixed 
ring  and  three  planets. 
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For  the  degenerate  natural  frequencies  (03=...  =6)n  i,  the  single-mode  instability 
boundaries  are  the  same  as  (5.28)  with 

r">  =  (/>“)’  =^[(k,psmlnc,f  +  (*„gsin lncr?  + 

2kspkipgsin(lm:s)sm(l7rc,.)cosln:(cs -cr+2y„)] 

Equation  (5.35)  applies  to  any  number  of  planets  N. 

For  the  combination  instability  of  degenerate  0)3=...  =C0n~i  and  distinct  ft).,  the 
combination  instabilities  (1=1)  of  a  distinct  mode  and  a  degenerate  mode  ahvays  vanish 
for  any  N  when  the  planet  meshes  are  in-phase.  A  more  general  condition  for  vanishing 
of  these  combination  instabilities  is  that  each  of  ( hs±l)/N ,  (lzr±l)/N  #  integer. 

As  an  example,  Figure  5.11  shows  the  instability  boundaries  for  a  planetary  gear 
with  three  equally  spaced  planets.  The  parameters  and  nominal  natural  frequencies  are 
given  in  Table  5.2  and  the  vibration  modes  are  shown  in  Figure  5.12.  The  mesh  phasing 
is  Ysn~Yni=0,  Y*r=l/2  and  the  contact  ratios  are  cs=1.4,  cr=1.6.  We  specify  £/=£?=£  (that 
is,  g=l ).  The  natural  frequencies  change  as  £  varies  because  the  average  mesh  stiffnesses 
depend  on  the  amplitude  of  the  mesh  stiffness  variation.  The  analytical  solutions  (solid 
lines)  from  (5.26),  (5.27),  (5.28),  and  (5.29)  agree  well  with  the  numerical  solutions  using 
Floquet  theory  and  numerical  integration.  Note  the  combination  instabilities  at  (O2+CO3, 
(O3+CO5  vanish  because  they  involve  the  distinct  ( Oh,  ay)  and  degenerate  ((03=0)4)  natural 
frequencies.  Figure  5.13a  shows  the  primary  instability  regions  for  different  contact 
ratios  cs,  cr;  the  stiffness  variation  amplitude  £=0.3  and  phasing  All  instabilities 

vanish  when  the  contact  ratios  are  integers.  The  size  of  the  2ah  instability  region  is 
primarily  affected  by  cr  while  insensitive  to  changes  in  cs.  This  is  because  the  dominant 
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deformation  of  mode  0?  occurs  in  the  ring-planet  meshes  (Figure  5.12).  The  264- 
instability  region  is  mostly  affected  by  cs  because  most  deformation  occurs  in  the  sun- 
planet  meshes  for  this  mode.  The  primary  instability  of  the  degenerate  modes  at  2(03=2(04 
is  influenced  by  both  cs  and  cr.  The  primary  instability  regions  are  maximized  for 
cs-cr=1.5,  which  is  expected  from  (5.33)  and  (5.35).  The  relative  phasing  ysr  between  the 
sun-planet  and  ring-planet  meshes  can  have  a  major  impact.  When  the  phasing  ysr=l/2 
and  other  parameters  are  the  same  as  in  Figure  5.13a,  the  2 0)3=2 0)4  instability  vanishes 
for  any  cs=cr  (Figure  5.13b).  This  is  because  T(I)  =0  in  (5.35)  for  ksp=krp,  g=l  and  cs=cr. 

The  above  analyses  show  that  contact  ratios  and  mesh  phasing  significantly  affect 
the  operating  condition  instability  regions.  In  practice,  particular  instabilities  can  be 
minimized  by  proper  selection  of  contact  ratios  and  mesh  phasing,  which  are  adjusted  by 
center  distance,  diametral  pitch,  pressure  angle,  tooth  addendum,  numbers  of  teeth,  and 
other  parameters. 


(2)  Sequentially  Phased  Planet  Meshes 

Here  we  consider  equally  spaced  planet  systems  where  the  sun-planet  and  ring- 
planet  meshes  are  sequentially  phased  with  ys„=(n-l)zs  /N  and  ym=(n-l)zr  /N.  This  case 
corresponds  to  zs  /N,  zr  /N#  integer  but  (zs+zr)/N= integer.  A  constant  phasing  ysr  exists 
between  the  sun-planet  and  ring-planet  meshes  for  each  planet.  For  the  sequential 
phasing,  the  Fourier  coefficients  in  (5.21)  satisfy 

- z: = 0  (5-36) 

when  lzs/N,  lzr/N  are  non-integer. 
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For  distinct  natural  frequencies  c$,  the  primary  instabilities  (1=1)  of  distinct  modes 
always  vanish  when  the  planets  are  sequentially  phased.  Physically,  it  means  that  these 
instabilities  are  not  excited  because  the  resultant  modal  excitations  from  the  sun-planet 
and  ring-planet  meshes  each  are  zero.  The  secondary  instabilities  (1=2)  of  distinct  modes 
also  vanish  when  2zs/N,  2zr/N  are  non-integer.  Similarly,  combination  instabilities  (1=1) 
of  distinct  (Op  and  ooq,  =  0  always  vanish  in  this  case. 

For  instabilities  involving  the  degenerate  modes,  it  is  difficult  to  obtain  simple 
expressions  for  the  instability  boundaries;  their  instability  conditions  can  be  calculated 
from  (5.28),  (5.30),  and  (5.31).  A  special  case  is  the  combination  instability  (1=1)  of 
degenerate  a>s=.. .=<%+/  and  distinct  ay  when  (zs±l)/N,  (zr±l)/N  each  is  non-integer.  In 
this  case,  the  combination  instabilities  of  distinct  and  degenerate  modes  vanish. 

Figure  5.14  shows  the  instability  boundaries  for  the  same  system  as  in  Figure  5.1 1, 
except  the  three  planets  are  sequentially  phased  with  %„=[0,  1/S,  2/3],  yn,=[0,  2/3,  1/3). 
Note  the  primary,  secondary,  and  combination  instability  regions  vanish  for  the  distinct 
natural  frequencies  Oh,  a%;  only  instabilities  involving  at  least  one  of  the  degenerate 
modes  003=0)4  exist. 

The  foregoing  results  for  equally  spaced  planets  are  summarized  in  Table  5.3. 


NASA/CR— 200 1  -210939 


92 


Figure  5.11  Instability  regions  for  the  system  in  Table  5.2  and  in-phase  meshes  Ysn=Y™ 
Ysr=l/2.  cs-1.4,  cr=1.6,  e=£i~£2  (g=l).  —  analytical  solution;  ***  numerical  solution. 


Figure  5.12  Mode  shapes  of  the  system  in  Table  5.2.  The  carrier  motion  is  not  shown. 
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Figure  5.13  Instability  regions  of  the  system  in  Table  5.2  for  different  contact  ratios. 
£=£i=£2=0.3  ( g-1 ).  The  planets  are  equally  spaced  with  in-phase  meshes  (ys„=y„,=0). 
(a)  ysr=0,  (b)  ysr=l/2. 
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Figure  5.14  Instability  regions  for  the  three-planet  system  in  Table  5.2  and  sequentially 
phased  planet  meshes  with  ys„=[0,  1/3,  2/3],  ym=[0,  2/3, 1/3].  cs=1.4,  cr=1.6,  £=£i=£2 
(g=l),  ysr=l/2.  —  analytical  solution;  ***  numerical  solution. 


Planet  Mesh 
Phasing 

Single-Mode  Instabilities 

Combination  Instabilities 

Distinct  Mode 

Degenerate 

Mode 

Distinct  + 
Distinct  Mode 

Degenerate  + 
Distinct  Mode 

In-phase 

from  (5.28), 
(5.30) 

=  integer 

N  N 

from  (5.26) 

from  (5.27) 

always  vanish 

Sequentially  phased 

— ^  integer 

N  N 

primary  always  vanish; 
secondary  from  (5.26) 
and  vanish  if  2zs  /N} 

2zr/N  ^  integer 

from  (5.28), 
(5.30) 

always  vanish 

from  (5.31)  and 
vanish  if  (zs±l)/N} 
(Zy±l )/N  &  integer 

Table  5.3  Instability  boundary  solutions  when  the  planets  are  equally  spaced  (satisfying 
(zs+zr)/N  =  integer). 
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2.  Unequally  spaced  planets 

When  the  planets  are  arbitrarily  spaced,  the  structured  modal  properties  do  not  exist 
and  the  general  expressions  for  instability  boundaries  cannot  be  further  simplified.  For 
the  practically  important  case  of  diametrically  opposed  planets,  however,  the  vibration 
modes  retain  the  well-defined  properties  (5.24)  and  (5.25).  The  following  discussion 
focuses  on  this  case  for  N/2  pairs  of  diametrically  opposed  planets.  For  the  sun-planet 
meshes,  each  pair  of  diametrically  opposed  planets  are  in-phase  ( yS(n+N/2 )  =Ysn)  for  even  zs 
and  are  counter-phased  (yS(n+N/2)  =T&»+l/2)  for  odd  zs.  Analogous  rule  applies  for  the  ring- 
planet  mesh  phasing.  Note  that  adjacent  planets  have  arbitrary  mesh  phasing  %„=y/,£r 
/(2n),  Ym=W^r/(2n),  n=l,...,N/2. 

When  the  sun-planet  and  ring-planet  meshes  are  both  counter-phased  (odd  zs,  zr), 
equation  (5.36)  holds  for  odd  /  because  a(JJ  =  -ct(s[l+NI2)  and  similar  relations  for 

aw>  b(s'J,  b^ .  Recalling  modal  property  (5.24),  =  E1^  =0  in  (5.32)  for  distinct  cq„ 

Cty  and  odd  /.  From  (5.26)  and  (5.27),  primary  and  combination  (1=1)  instabilities  of 
distinct  modes  always  vanish  in  the  counter-phased  case.  Physically,  these  instabilities 
are  eliminated  because  the  modal  excitations  from  each  pair  of  diametrically  opposed 
planets  always  cancel  each  other.  The  instability  regions  involving  degenerate  modes  are 
obtained  from  numerical  evaluation  of  the  eigenvalues  of  (5.30). 

When  pairs  of  opposing  sun-planet  or  ring-planet  meshes  are  in-phase  (even  zs  or 
even  zr),  no  simple  expressions  for  instabilities  regions  are  available;  numerical  solutions 
are  obtained  from  (5.26),  (5.27),  (5.30),  and  (5.31). 
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5.2.4  Dynamic  Response  and  Contact  Loss 

When  planetary  gears  are  operated  inside  an  instability  region,  damping  and 
nonlinearities  from  friction,  tooth  separation,  etc  bound  the  unstable  linear  model  motion. 
Figure  5.15a  shows  the  RMS  steady-state  planet  response  amplitude  versus  mesh 
frequency  for  the  same  system  as  in  Figure  5.11  and  the  stiffness  variation  £=0.3. 
Rayleigh  damping  C=(0.07)*M+(0.07)*K  is  added  to  system  (5.2)  and  the  force  vector 
is  ¥=[-2000  1000  0  0  0]T N.  The  solutions  are  obtained  from  numerical  integration  using 
mesh  stiffnesses  in  rectangular  waveforms  (Figure  5.15b).  The  degenerate  modes  (03-0)4 
only  have  planet  motion,  so  they  are  not  excited  for  the  torques  applied  to  the  carrier  and 
sun.  The  combination  instability  (02+0)5  is  sufficiently  damped  so  that  it  is  not  apparent 
in  Figure  5.15a  and  16a.  The  resonance  excited  by  the  primary  instability  2 0)5  is 
extremely  large  because  tooth  separation  is  not  considered;  the  mesh  stiffnesses  are  pre¬ 
specified  functions  of  time  (Figure  5.15b).  In  practice,  tooth  separation  (clearance 
nonlinearity)  occurs  for  large  dynamic  responses  and  its  effects  are  dramatic.  Figure 
5.16a  shows  the  response  for  the  same  system  as  in  Figure  5.15a,  but  tooth  separations  is 
modeled.  The  mesh  stiffness  ks„  or  km  is  set  to  zero  if  the  corresponding  tooth 
deformation  Ssn<0  or  Sm<0  at  any  step  of  the  integration.  The  response  amplitude  of  the 
2(05  primary  instability  is  significantly  reduced  from  that  in  Figure  5.15a.  Moreover,  a 
softening  jump  phenomenon  occurs.  Sun-planet  tooth  separation  (ksp=0)  is  apparent  in 
Figure  5.16b  for  Q~2a £.  The  interactions  of  mesh  stiffness  variation  and  clearance 
nonlinearity  was  studied  by  Kahraman  and  Blankenship  (1996,  1997)  for  single-mesh 
gears.  Their  effects  on  multi-mesh  planetary  gears  need  additional  investigation. 
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Figure  5.16  (a)  Steady  state  RMS  of  dynamic  planet  rotation  versus  mesh  frequency  £2 
when  tooth  separation  is  considered.  The  parameters  are  as  in  Figure  5.11  with  £=0.3. 
Circles  (o)  indicate  increasing  speed  £2  and  crosses  (x)  indicate  decreasing  £2.  (b)  Sun- 
planet  tooth  separation  ( ksp=0 )  occurs  for  £2=4.5kHz~2(Q$. 
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CHAPTER  6 


SUMMARY  AND  FUTURE  WORK 


6.1  Summary  and  Benefits 

This  work  analytically  investigates  several  key  issues  in  planetary  gear  dynamics. 
The  main  results  and  their  benefits  are  summarized  for  each  specific  topic. 

1.  Dynamic  Model  of  Planetary  Gears 

A  lumped-parameter  model  is  developed  for  spur  planetary  gears.  Critical  factors  to 
gear  vibration  are  considered,  including  arbitrary  planet  spacing  and  phasing,  gyroscopic 
effects  that  are  important  in  high-speed  applications,  mesh  stiffness  variation,  and 
transmission  error  excitation.  The  model  is  applicable  to  general  epicyclic  gears  with 
various  configurations.  It  is  suitable  for  the  dynamic  analysis  of  critical  issues  in 
planetary  gear  vibration  and  useful  for  design  guidance  (response,  natural  frequencies, 
etc.).  This  model  is  a  key  building  block  that  can  be  expanded  to  couple  with  the  housing 
and  include  multiple  stages. 

2.  Modal  Properties 
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This  study  characterizes  the  natural  frequency  spectrum  and  vibration  mode 
properties  in  planetary  gears.  Planetary  gears  possess  rich  modal  structure  that  is  crucial 
for  eigensensitivity,  dynamic  response,  and  stability  analysis.  Rigorous  characterization 
of  these  special  properties  is  a  fundamental  advance  for  planetary  gear  vibration  research. 
These  properties  are  capsulized  as  follows: 

•  When  N  planets  are  equally  spaced,  the  vibration  modes  can  be  classified  into  six 
rotational  modes  with  distinct  natural  frequencies,  six  pairs  of  translational  modes 
with  degenerate  natural  frequencies  of  multiplicity  two,  and  three  groups  of  planet 
modes  with  degenerate  natural  frequencies  of  multiplicity  N-3.  Each  type  of  vibration 
mode  has  unique  properties  due  to  the  cyclic  symmetry  of  the  system. 

•  When  N/2  pairs  of  planets  are  diametrically  opposed,  the  rotational  and  planet  modes 
have  the  same  structure  as  for  equally  spaced  planets;  translational  modes  lose  their 
degeneracy  but  retain  their  distinct  properties. 

•  When  the  planets  are  arbitrarily  spaced,  the  rotational  and  translational  modes 
generally  lose  their  well-defined  structure;  the  remarkable  properties  of  planet  modes 
are  not  affected  by  planet  spacing. 

3.  Design  Parameter  Variations 

The  effects  of  design  parameter  variations  on  planetary  gear  free  vibration  are 
analytically  investigated.  The  sensitivities  of  natural  frequencies  and  vibration  modes  to 
key  parameters  are  expressed  in  simple,  closed-form  formulae  that  allow  one  to  quickly 
identify  the  parameters  that  most  impact  the  modal  properties.  Well-defined  veering  rules 
are  derived  to  examine  the  dramatic  changes  of  vibration  modes  that  are  possible  even  for 
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small  parameter  variations.  These  results  rely  on  the  structured  vibration  mode  properties 
characterized  before.  For  design  use,  the  results  allow  qualitative  and  quantitative 
assessment  of  the  effects  of  design  parameters  on  planetary  gear  free  vibration: 

•  Each  of  the  sun-planet  and  ring-planet  mesh  stiffnesses  affect  one  rotational  mode, 
one  translational  mode,  and  one  planet  mode. 

•  The  transverse  support  stiffnesses  and  component  masses  of  the  carrier,  ring,  and  sun 
only  affect  translational  modes. 

•  The  torsional  support  stiffnesses  and  moments  of  inertia  of  the  carrier,  ring,  and  sun 
only  affect  rotational  modes. 

•  The  planet  bearing  stiffness  and  inertia  are  critical  parameters  that  affect  most  modes. 

•  The  operating  speed  does  not  affect  rotational  and  planet  modes;  translational  mode 
natural  frequencies  at  zero  speed  split  into  distinct  ones  as  speed  increases. 

4.  Mesh  Stiffness  Variation  Excitation 

The  parametric  instabilities  excited  by  mesh  stiffness  variations  are  investigated  for 
two-stage  and  planetary  gear  trains.  The  operating  conditions  leading  to  parametric 
instabilities  are  analytically  determined  and  numerically  verified.  The  following  findings 
are  obtained  from  this  study. 

•  The  contact  ratios  and  mesh  phasing  strongly  impact  the  instabilities  induced  by  mesh 
stiffness  variation.  Simple,  exact  formulae  derived  from  perturbation  analyses  provide 
design  guidance  to  suppress  particular  instabilities  by  adjusting  the  contact  ratios  and 
mesh  phasing. 
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•  For  two-stage  gears  in  a  countershaft  configuration,  the  mesh  excitations  interact 
when  one  mesh  frequency  is  an  integer  multiple  of  the  other  and  dramatically  change 
the  instability  conditions  compared  to  two  decoupled  excitations. 

•  For  planetary  gears,  the  structured  modal  properties  lead  to  vanishing  of  certain 
instabilities,  including  the  combination  instability  of  distinct  and  degenerate  natural 
frequencies  when  planet  meshes  are  in-phase,  and  the  single-mode  and  combination 
instabilities  of  distinct  natural  frequencies  when  planet  meshes  are  sequentially 
phased.  Along  with  contact  ratios,  selection  of  planet  phasing  is  an  effective  means  to 
avoid  parametric  instability  and  its  associated  large  response. 

6.2  Future  Work 

To  thoroughly  understand  planetary  gear  dynamics  and  develop  reliable  design 
tools  for  noise  and  vibration  reduction,  this  research  needs  to  be  advanced  in  several 
challenging  areas  discussed  below. 

1.  Nonlinear  Effects 

Tooth  separation  nonlinearity  strongly  affects  the  dynamic  behaviors  of  single-mesh 
gears,  including  jump  phenomena  and  multiple  steady  state  solutions  (Blankenship  and 
Kahraman  1995;  Kahraman,  1992;  Kahraman  and  Blankenship,  1996,  1997;  Kahraman 
and  Singh,  1991;  Parker  et  al.,  2000b;  Rook  and  Singh,  1995;  Theodossiades  and 
Natsiavas,  2000).  Typical  methodologies  such  as  direct  integration,  harmonic  balance, 
perturbation,  and  shooting  method  were  used  in  the  above  analyses.  This  nonlinearity  has 
not  been  fully  investigated  in  planetary  gears.  The  contact  loss  nonlinearity,  coupled  with 
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mesh  stiffness  variation,  transmission  error  excitation,  and  mesh  phasing,  is  expected  to 
impact  dynamic  responses,  load  sharing  among  planets,  bearing  loads,  tooth  fatigue,  and 
gear  noise.  Analytical  study  of  nonlinear  effects  will  provide  improved  design  guidance 
for  planetary  gears  based  on  dynamic  response. 

2.  Transmission  Error  Modeling 

For  single-mesh  gears,  static  transmission  error  (STE)  is  naturally  defined  and 
widely  used  as  a  metric  for  vibration  and  noise  (Gregory  et  al.,  1963-64;  Smith,  1987; 
Ozguven  and  Houser,  1988b).  The  specification  of  transmission  errors  in  planetary  gears 
is  complicated  due  to  the  carrier  rotation  and  multiple  meshes  coupled  together.  Donley 
and  Steyer  (1992)  simply  used  a  “net  transmission  error”  that  is  the  average  of  all 
individual  transmission  errors  of  the  sun-planet  and  ring-planet  meshes.  Although  this 
simplification  reduces  the  computational  and  modeling  effort,  it  is  a  significant 
simplification  difficult  to  interpret  of  the  dynamic  mesh  forces.  Other  studies  (Kahraman, 
1994b,  1999;  Kahraman  and  Blankenship,  1994)  used  Fourier  series  to  represent  STE  in 
each  mesh,  but  did  not  address  the  physical  modeling  related  to  profile  error,  pitch  error, 
run-out  error,  and  misalignment.  It  remains  a  question  how  to  synthesize  transmission 
errors  from  these  parameters  in  planetary  gears.  More  investigations  are  needed  to  define, 
calculate,  and  measure  a  useful  analog  of  STE  in  planetary  gears. 

3.  Dynamic  Response 

Dynamic  response  from  analytical  models  is  useful  to  predict  the  noise  and 
vibration  early  in  the  design  stage.  To  reliably  determine  dynamic  response,  it  requires 
precise  representation  of  dynamic  excitation  from  tooth  meshes.  The  time-varying  mesh 
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stiffness  and  transmission  error  for  single-mesh  gears  can  be  calculated  from  the  contact 
analysis  program  LDP  (Houser,  1990)  or  other  tools.  A  major  challenge  is  to  match 
dynamic  response  obtained  from  analytical  models  using  these  approximations  as  input 
with  benchmark  computational  and  experimental  results.  No  such  comparison  exists  for 
planetary  gears,  and  this  restricts  the  application  of  analytical  models  in  practical  design. 
In  addition,  the  effects  of  major  design  factors  ( e.g .,  contact  ratio,  mesh  phasing,  mesh 
stiffness  variation,  support/bearing  stiffness,  tooth  modification)  on  dynamic  response 
need  be  characterized  for  design  guidance.  This  investigation  will  start  with  two-stage 
spur  gears  and  then  extend  to  more  complicated  planetary  gears. 

4.  Ring  Gear  Flexibility 

The  ring  body  is  considered  rigid  in  this  dissertation.  In  practice,  internal  gears  with 
thin  rims  may  distort  elastically.  Experimental  measurements  (Ma  and  Botman,  1984) 
and  finite  element  computations  (Kahraman  and  Vijayakar,  2000)  have  shown  dramatic 
changes  in  load  sharing  among  planets  from  the  ring  flexibility.  Because  the  ring  is 
usually  connected  to  the  housing,  its  flexibility  directly  influences  the  transmission  of 
structure-borne  noise.  The  elastic  distortion  of  the  ring  can  also  affect  the  mesh 
stiffnesses  and  lead  to  contact  loss.  Qualitative  and  quantitative  analyses  of  the  effects  of 
ring  flexibility  is  necessary  to  improve  understanding  of  physical  behaviors.  The 
fundamental  task  is  to  extend  the  current  discrete  model  by  including  the  continuous 
deformation  of  the  ring  body.  The  analytical  and  computational  methods  discussed  in 
this  work  can  likely  be  used  on  this  extended  model. 

5.  Experimental  Verification 
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5.  Experimental  Verification 

The  inherent  complexity  of  gear  dynamics  requires  experimental  validation  of 
analytical  and  numerical  results.  For  single-mesh  gear  models,  there  are  extensive 
experiments  verifying  analytical  findings  as  well  as  presenting  questions  for  research  and 
design.  The  analytical  models  of  planetary  gears,  however,  have  not  been  properly 
validated  because  of  the  lack  of  benchmark  experimentation.  The  calculated  vibrations 
from  the  finite  element  tool  Calyx  agree  well  with  measured  results  at  NASA  Glenn 
Research  Center  (Krantz,  1992)  and  General  Motors  (Blankenship  and  Kahraman,  1996; 
Kahraman  and  Blankenship,  1997,  1999a,  b;  Parker  et  al.,  2000b).  More  experimental 
studies  are  needed  to  develop  reliable  dynamic  models  suitable  for  practical  design  and 
investigate  critical  issues  in  planetary  gear  vibration. 
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APPENDIX  A:  SYSTEM  MATRICES 


M  =  diag(Mc 

Mft  =  diag(mh ,mhJh  / r2 ),  h  =  c,r,s,l,-",N 
0  -  2  mh  0 
G  =  2  mh  0  0 

0  0  0 

Kn  =  diag{mc ,  mc  ,0,  mr,  mr  ,0  ,ms,  ms  ,0,  mp,mp,0,---,mp,  mp  ,0) 
Kj  =  J/ag(Kci,Kri,Kj6,0,---,0) 

Khb  =  diag(khx ,  khy ,  khu ),  h  =  c,r,s 

ZK"i  0  0  K*2  K2c2  -  Knc2 

Zk”i  0  k'2  K^2  •••  Kf2 

Sjy'W  IS  ^ 

*^sl  *^s2  *^s2  ”•  ^2 

“m 

KpP  0  -  0 

symmetric  : 


1  0  -  sin  y/n 

-cos  y/n 

sin  Wn 

0" 

w 

^  a 

11 

1  cos  wn  K  nc2  = 

-sin  Wn 

-cos  y/n 

0 

symmetric  1 

0 

-1 

0 

K”3  =  diag(kpn  kpn  0) 
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k;.'i = mo 


sin  y/n,  -  cos  y/n,  sin  y/ni 


cos  y/„ 


symmetric 


-sin  Vm 
cos  y/n, 

1 


-  sin  ^rasinorr 


Ki=kM 


cos^n7sinar 


sin  aT 


-  sin  y/n1  cos  or,.  sin  y/rn 

-  cos  ^,77cos  aT  -  cos  y/rn 

-cosaT  1 


Kr"3  =*„,(') 


sin2  cr. 


-  cos  a,,  sin  ar 
cos  2ar 


symmetric 


-  sin  ar 
cos  ar 
1 


k;,  =  **(o 


sln“  Wrn  -  C0S¥sn  SU1  Wsi 


cos  y/sl 


symmetric 


- sin  ¥s„ 

cos  ys„ 

1 


K  ns2  =  ks„(t) 


sin  ^sina,, 
-cos^„sinas 


sin  Wsn  cos  as 
-cosyrs„cosas 


-  siners  -  cos  as 


■  sin  Wsn 
cosy/sn 
1 


k?3=MO 


sin2  a * 


cosas  sin  as 
cos  2as 


symmetric 


-sinc^ 
-  cos  as 
1 


T(/)  =  [0,0, 7;  ,0.0. 7;,o,o,7;,o,-o]r 

F(0  =  [0,Fr,FJ,F1  •••Fw]r, 

F,  =  [sin  y/m  -  cos  ysm,l]T , 

F,  =  k,nesn  [sin  y/sn  -  cos  y/sn  ,\]T 

F„  =^„,[sinar,-sin«r,-lf  +^„eJ„[sina,,-sin«r,l]r 


Note  that  the  mesh  stiffness  variation  is  modeled  by  the  time-varying  stiffnesses 
ksn(0  and  km(t). 
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APPENDIX  B:  EXAMPLE  SYSTEM 


The  planetary  gear  used  in  the  U.S.  Army’s  Helicopter  OH-58  is  the  benchmark 
example  in  the  work.  Two  cases  of  the  system  are  used  in  the  study. 

Case  I:  The  ring  gear  is  free  and  the  system  has  3(N+3)  degree  of  freedom,  where  N  is 
the  number  of  planets.  The  model  parameters  are  given  in  Table  C.l. 

Case  II:  The  ring  gear  is  fixed  (kr=  kru=  0)  and  other  parameters  are  the  same  as  Table 
C.l.  The  system  has  four  planets  and  18  degrees  of  freedom. 


Sun 

Ring 

Carrier 

Planet 

Mass  (kg) 

0.4 

2.35 

5.43 

0.66 

I/r2  (kg) 

0.39 

3.00 

6.29 

0.61 

Base  diameter  (mm) 

77.42 

275.03 

177.8 

100.35 

Teeth  number 

27 

99 

35 

Mesh  stiffness  (N/m) 

ksp  =  kip—  km=  5  x10s 

Bearing  stiffness  (N/m) 

°o 

II 

II 

II 

A? 

II 

Torsional  stiffness  (N/m) 

Ku  =  109  ksu=  kcu  =  0 

Pressure  angle  (°) 

Os  =  0?  =  a  =  24.6 

Table  B.lModel  parameters  of  the  planetary  gear  in  the  U.S.  Army’s  helicopter  OH-58 
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